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Abstract 

The semiclassical (zero-dispersion) limit of the one-dimensional focusing Nonlinear Schrodinger equation 
(NLS) with decaying potentials is studied in a full scaling neighborhood D of the point of gradient 
catastrophe {xo,to). This neighborhood contains the region of modulated plane wave (with rapid phase 
oscillations), as well as the region of fast amplitude oscillations (spikes). In this paper we establish the 
following universal behaviors of the NLS solutions near the point of gradient catastrophe: i) each spike has 
the height 3\qQ{xQ, to, e) | and uniform shape of the rational breather solution to the NLS, scaled to the size 
0(e); ii) the location of the spikes are determined by the poles of the tritronquee solution of the Painleve 
I (PI) equation through an explicit diffeomorphism between D and a region into the Painleve plane; iii) 
if {x,t) e D but lies away from the spikes, the asymptotics of the NLS solution q{x,t,e) is given by the 
plane wave approximation qo{x, t, e), with the correction term being expressed in terms of the tritronquee 
solution of PI. The latter result confirms the conjecture of Dubrovin, Grava and Klein [12] about the 
form of the leading order correction in terms of the tritronquee solution in the non-oscillatory region 
around (a;o,to)- We conjecture that the PI hierarchy occurs at higher degenerate catastrophe points and 
that the amplitudes of the spikes are odd multiples of the amplitude at the corresponding catastrophe 
point. Our technique is based on the nonlinear steepest descent method for matrix Riemann-Hilbert 
problems and discrete Schlesinger isomonodromic transformations. 
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1 Introduction and main results 



In this paper we consider the focusing Nonhnear Schrodinger (NLS) equation 

iedtq + e^dlq + 2\q\^q = 0,^ (1-1) 

where x S M and t > are space-time variable and e > 0. It is a basic model for self-focusing and 
self- modulation, for example, it governs nonlinear transmission in optical fibers; it can also be derived as 
a modulation equation for general nonlinear systems. It was first integrated by Zakharov and Shabat [30] 
who produced a Lax pair for it and used the inverse scattering procedure to describe general decaying 
solutions (lim|^|_j.oo q{x, 0) = 0) in terms of radiation and solitons. Throughout this work, we will use the 
abbreviation NLS to mean "focusing Nonlinear Schrodinger equation" . 

Our interest in the semiclassical (zero-dispersion) limit (s — > 0) of NLS stems largely from its mod- 
ulationally unstable behavior. As shown by Forest and Lee [14], the modulation system for NLS can be 
expressed as a set of nonlinear PDE with complex characteristics; thus, the system is ill posed as an 
initial value problem with Cauchy data. As a result, initial data 

qix, 0, e) = A(x)e^*(^)/^ $(0) = 0, (1-2) 

i.e. a plane wave with amplitude modulated by A{x) and phase modulated by ^(x) are expected to 
break immediately into some other, presumably disordered, wave form when the functions A{x) and ^(x) 
possess no special properties. 

However, in the case of an analytic initial data, the NLS evolution displays some ordered structure 
instead of the disorder suggested by the modulational instability, as shown by Fig. 1 and Fig. 2, see 
[23], [28] and [8]. These figures clearly identify the regions where different types of behavior of the 
solution q{x, t, e) to the NLS appear. These regions are separated by some curves in the x, t plane that 
are asymptotically independent of e. They are called breaking curves or nonlinear caustics. Within 
each region, the strong asymptotics of q{x,t,e) can be expressed in terms of Riemann Theta-functions 
to within an error term of order 0{e) that is uniform on compact subsets of the region ([19] for the pure 
soliton case and [26] for the general case). In this context, regions of different asymptotic behavior of 
q{x,t,e) corresponds to the different genera of the hyperelliptic Riemann surface TZ{x,t), whose Theta- 
functions enter in the asymptotic description. Note that Ti{x, t) does not depend on e. In the very 
first (genus zero) region (that contains t — 0), the asymptotics of q{x,t,e) is expressed through the 
branch-point a{x,t) oiTZ{x,t) (the branch-points are Schwarz-symmetrical) as 

qo{x,t,e) = A(x,Oei*(^'*) = c^^^^^ ^^^-J /<-'/{25Ja(C,t)de+[4(5Ra(.,.))^-2(c,a(.,r))^]d.} ^ 4 (^.g) 

^Note that this equation differs by the coefficient 2 from the NLS equation considered in previous papers [26]-[27], whose 
results we use here; this is the correct equation for the particular form of the time evolution of the scattering data (see 
(2-5)), adopted in [26]-[27] and in the present paper. Alternatively, one can keep the "old" form of the NLS, but scale by 2 
the time variable while using the results of [26]-[27]. 

^This formula was proven in [26], but, as stated in Theorem 1.1 there (with only part of the differential form), is 
correct only for t = 0. 
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Figure 1: Absolute value \q{x, t, e) \ of a solution q{x, t, e) to the focusing NLS (1-1) versus x, t coordinates. 
Here A{x) = e^^ , = tanhx and e = 0.03. 

where the error estimate is uniform on compact subsets of the genus zero region (away from the breaking 
curve) . 

A detailed study of asymptotic behavior of q{x, t, e) along the first breaking curve (a neighborhood of 
the tip of this curve was excepted) together with error estimates were conducted in our previous work [4] . 
Considering one of the pieces of the breaking curve (to the left, a; < 0, or to the right, a; > 0, from the tip, 
see Fig. 2), we introduced scaled coordinates S{x,t) = ^'d{x,t) + ^K{x,t)e\ne, that map a neighborhood 
of the left (right) breaking curve onto a horizontal strip in C, containing negative (positive) real semi-axis. 
In the coordinates z9, k, the oscillatory region of the left breaking curve is in the upper 5'-half-plane, and 
the shape of oscillations in the S-plane is shown on Fig. 3. 

The tip-point XQ,to of the breaking curve is called a point of gradient catastrophe, or elliptic 
umbilical singularity ([12]).) It is evident from the numerical simulations shown on Fig. 1 and Fig. 
2, as well as from the analysis of the spectral plane (presented below), that the behavior of the solution 
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Figure 2: A{x) = e"^ , <i>'(x) = tanha; and e = 0.03. 



q{x, t, e) at the tip of the breaking curve is very different from the behavior elsewhere on the breaking 
curve. 

The main goal of this paper is to analyze the leading order asymptotic behavior of the solution 
q{x, t, e) on and around this special point of transition. More precisely we will examine a neighborhood 
D of (xq, to) that is shrinking at the rate 0{e^) as e — > 0. As the first step, we construct a diffeomorphism 
V = v{x, t, e), that maps D onto a bounded disk y e C, where V is independent of e and v{xq, to, e) = 0. 
It turns out that the leading order behavior of q(x, t, e) in D can be conveniently described through a 
specific tritronquee solution (see Section 4.2.1) to the Painleve I (PI) equation 

y" = - V. (1-4) 

That is why the diffeomorphism v(x,t,e) and the complex w-plane will be referred to as the Painleve 
coordinatization of D and the Painleve plane respectively. Note that the diffeomorphism v{x,t,e) near 
the point of gradient catastrophe plays similar role to the map S{x, t, e) for the rest of the breaking curve. 

The point of gradient catastrophe (xo,to) can be defined in terms of the physical variables {x,t) as 
the point, where the genus zero approximation (1-3) of the solution q{x,t,e) to (1-1), (1-2) develops an 
infinite x-derivative in A{x,t) and/or in ^x{x,t) (while A{x,t) and ^x{x,t) stay finite). In terms of the 
time evolution of the spectral data (see Fig. 4 below), the point of gradient catastrophe is the point 
where the birth/collapse of the new main arc (band) happens exactly at the end of an existing main arc 
(see Fig. 4). 

As it was mentioned above, the asymptotics of q{x,t,e) within the region of a given genus is given 
explicitly in terms of the Ricmann Theta-functions, associated with TZ{x,t), with the accuracy 0{e), see 
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Figure 3: The graphs depicts the ration ^^^^^^1-^ in a vicinity of a breaking curve point a: = 0.1 , t = 
0.1570482549 with e = 10"^". Here go is the genus zero approximate solution, q is the leading order 
approximation near the breaking curve, constructed in [4]. The graph is shown in the variables 9,x 
and the scale is not uniform: x measures distances in the scale e|lne| whereas 9 in the scale e. The 
distance between consecutive "ranges" , should be much longer than the separation of the neighboring 
peaks within the same range. The typical size of the hills is e in all directions while the separation in the 
longitudinal direction is e\ lne|. This picture is consistent with a full-blown genus-two regime, where the 
solution is quasiperiodic at the scale e: as we progress into the genus-2 region the separation reduces to 
the natural scale e. If we were to plot this in the (x, t) plane the only difference would be a linear change 
of coordinates and the "ranges" would be parallel to the breaking curve. 

[26], [19]. (In the case of genus zero, see (1-3).) The leading order behavior of q{x,t,e) in order O(elne) 
strips around left and right branches of the breaking curve, found in [4], has the accuracy O(y^). It 
was conjectured by B. Dubrovin, T. Grava and C. Klein in [12] that the error estimate near the gradient 
catastrophe is of order 0(e^/^) and is expressible in terms of the tritronquee solution of PI (see Remark 
5.2 for more details on their conjecture). 

1.1 Description of results 

The results of this paper not only confirm the conjecture of Dubrovin, but in fact go beyond it, as 
we provide the leading order behavior together with the accuracy estimate in the whole domain D 
around the point of gradient catastrophe xq, to, that include the oscillatory part of D. The following 
notations are useful in describing our results. 

Let Vp — {fp,i, • • • ,fp,Ar} denote the set of poles of the tritronquee solution y{v) in V , i.e, Vp C V. 
Let Bs_j denote the disk of radius S > centered at Vp j, j = 1, ■ ■ ■ , N, Bs = Bsj and Ks = V \ Bg. 
Denoting a + ib — a{xo, to), we prove the following facts: 



1. There is a one to one correspondence between the poles of the tritronquee solution y{v) 
within V and the spikes of the NLS solution q within D. Each spike is centered at the corresponding 
(xp^j,tp,j) = t;"i(wpj), where 

[2b / 2 \ 

v{x,t,e) = —^^\ — [x - xo + 2{2a + ib){t ~ to)] [1 + 0{e^)) (1-5) 

£5 V O \ / 

uniformly in D, with the nonzero constant C explicitly defined by (3-44) in terms of the scattering 
data; 

2. Each spike has the fixed height of 3|(7o(a^o, ^o)! + 0{s^^^), i.e., the height of each spike is three 
times the amplitude at the gradient catastrophe, see Theorems 6.2, 6.3; 

3. Each spike has the universal shape of the (scaled) rational breather solution to the NLS 

eq. (1-1), see Fig. 10, i.e, 

qix^e) = et*(--*^)Q,. ^^Iea^ (i + o{ei)), (1-6) 

where the rational breather 



\ l + 462(^ + 4ar/)2 + 16feV/ 



satisfies the NLS eq. (1-1) with space-time variables ^, rj. This breather approximation of the spike 
is valid in the domain Bj — v^^{Bs_j) of {xp,j,tpj), where S = 0{e^^^), see Theorems 6.2) and [4]. 
The size of each spike in the physical plane (the size of Bj) is thus 0(e), which is consistent with 
the size of spikes along the breaking curve (away from {xo,to)) and within the bulk of the genus 
two region, see above. Remarkably, the breather-like shape of spikes (see Fig. 3) is also preserved 
as one moves along the breaking curve away from the point of gradient catastrophe. The two zeroes 
("roots") and the maximum of each breather, shown on Fig. 10, occur at the same time (within 
the accuracy of our approximation). 

4. In Thm. 5.2 we show that if 5 > is a small fixed number then 

q{x,t,e) = (5-2et5(4r) + 0{ei)) x 



exp — 



i$(xo,io) - {aix - xo) - (2a2 - b^){t ~ to)) + etSR (y^i?/(t^))] (1-8) 



uniformly in Kg = v^^{Ks), where C is a nonzero constant explicitly given by (3-44), y{v) = 
y{v{x,t,e)) is the tritronquee solution and Hj — ^{y'{v))^ + vy{v) — 2y^{v). Equation (1-8) is 
consistent with the conjecture of [12] (see Remark 5.2). 
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5. If (5 = Ole"), where p e (0, g), and Ks = v ^{Kg), then equation (1-8) will be uniformly valid in 
Kg provided that in the error term will be replaced by ei~'^, i.e. 

qix,t,e) = (&-2et3(4^) +0(et-3'')) x 
i$(xo,io) - ia{x - xo) - (2a2 - b^){t - to)) + ei^ • (1-9) 



expf 



Note that - since y{v) has a double pole and Hiiv) a simple pole - the term y{v) is actually of 
order e~'^'^ and Hj{v) of order e~'^; clearly the description in terms of the tritronquee cannot be 
pushed "too close" to the pole/spike. 

We also make the Conjecture 6.1 that the amplitudes of the spikes near any gradient catastrophe 
point in the genus zero phase are odd multiples of the amplitude at the point itself. In addition we can 
speculate that the shape of the spikes at the higher-degeneracy cases should be related to the higher 
rational breather recently investigated in [1]. 

Among other results obtained in this paper we mention the proof that the two branches of the breaking 
curves form a corner (wedge) at the point of gradient catastrophe {xq, to) and give explicit expression, see 
(3-58), of the angle between the breaking curve in terms of C and a{xo,to)- We further prove that the 
diffeomorphism v{x,t,e) maps this corner in to the sector ^ < argu < ^ of the complex u-plane, see 
Fig. 9. This is consistent with another conjecture, stated in [12]: all the poles of the tritronquee solution 
y{v) are contained within the sector ^ < argu < which is a longstanding question in the theory of 
Painleve equations. According to our results, the set of spikes near the point of gradient catastrophe is, in 
fact, the visualization of the poles of the tritronquee solution y{v) to PI. For example, it was calculated 
in [18] that the very first pole Wp_i of y{v) has argtip^i = ^ and \vp^i\ ~ 2.38. Using (1-5), we calculate 
{xp,i,tp^i) — (0, 2.79126) for the NLS evolution of the initial data q{x, 0, e) = sechx. Numerical simulation 
of this evolution with e — shows the first spike at tg « 2.8, which is in a very good agreement with 
tp^i, see Example 6.1 and Fig. 11 

Our results are based on the nonlinear steepest descent analysis of the inverse scattering transform of 
the semiclassical NLS (1-1) that was set in the form of a matrix Riemann-Hilbert problem (RHP). This 
RHP depends on the semiclassical parameter e and external parameters x, t. As it is now well known, the 
nonlinear steepest descent method reduces this RHP into what is called a model RHP, which has some 
piece- wise constant jump matrices on a set of arcs in C, known as main arcs or bands. Solution to the 
model RHP represents the leading order term of the solution to the original RHP, and the main arcs are 
exactly the branch-cuts of the hyperelliptic Riemann surface TZ{x,t), associated with the problem. Note 
that, due to the symmetries of the focusing NLS, the set of main arcs is always Schwarz-symmetrical. 
The asymptotic behavior of the solution q{x,t,e) in a vicinity of {x,t) as e — > profoundly depends on 
the genus of TZ{x,t), as seen on Fig. 2 and Fig. 1, with the change of genus of Tl{x,t) corresponding 
to the transition between different asymptotic regimes of q{x,t,e). In these terms, the point of gradient 
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catastrophe can be described as the point where a new main arc (a new branch-cut oiTZ{x,t)), together 
with its complex-conjugated main arc, branches out from the endpoint of the existing main arc. 

While reducing the original RHP to the model RHP, the error is controlled through the use of local 
parametrices, which are designed to "mimic" the local behavior of the solution at the points of interest, 
for example, at the endpoints of main arcs. At the regular points {x,t) in the physical plane, these 
parametrices can be constructed through Airy functions. As we show in Section 4, at the point of gradient 
catastrophe (a^Oj^o): the parametrix V is constructed through the tritronquee solution y{v) of the PI. 
It is well known that Painleve equations can be expressed as conditions of isonionodromic deformations 
for certain rational 2x2 systems of ODEs with rational coefficients [17, 15]. The parametrix V at the 
point of gradient catastrophe is built through the fundamental solution to the system of ODEs 

associated with the PI. The occurrence of parametrices build out of Painleve associated linear systems 
is not unexpected here, as they often appear in various RHPs related to random matrices/orthogonal 
polynomials, for example: 

• The case of random matrices with a soft-edge where the density vanishes to order (z — q)^'^^-'^/^, 
corresponding to the even PI hierarchy [10]; 

• The case of random matrices where a spectral band splits into two (P2 equation [5] and hierarchy) ; 

• The trailing edge of the region of oscillations in the small-dispersion limit of KdV [9] (PH equation). 

The novelty of our work lies in the fact that the matrix v), and thus, the parametrix V, is not 
defined (has poles) at the poles Vp £ V of the tritronquee solution y{v). To our best knowledge, this 
paper contains the first example of parametrices with singularities, that were successfully used to control 
the errors. 

Solving the RHP at or near the pole Vp of y{v), i.e., studying the shape of the spike, require several 
additional steps, which can be briefly listed as: 

• Factorization ^'(^; v) — G'(^, u)*I'(^, u), where G(^, v) is a "simple" matrix with singularity at w = Vp, 
and \l/(^,i;) is regular aX v — Vp. This factorization was introduced by D. Masoero in [22]. The 
existence of the limit of ^(^, w) as w — > Vp that is uniform in a certain region of the spectral ^-plane 
(see Appendix A), is an important part in establishing the shape of spikes; 

• Construction of ^i, solution of a modified model RHP, that is needed at the poles v — Vp of the 
tritronquee solution y{v). was obtained by means of the discrete Schlesinger isomonodromic 
deformations of ^E'o. This type of deformations were used, for example, in [3, 4]; 

• Construction of the new parametrix Vi for the modified model RHP, using 5'(^,w). The formula 
36 -I- 0{e^^^) for the height of spikes follows immediately from Vi] 
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• Additional transformation of ^i, called partial Schlesinger transformation, is used to obtain the 
shape of the spikes, see (1-6). 

Finally, it is clear that our method can be modified to handle higher order (degenerate) gradient 
catastrophes, where k main arcs, fc > 1, simultaneously emerge from the endpoints of an existing main 
arc. The parametrices in these cases can be written in terms of the higher members of the PI hierarchy, 
and that one should expect the height of the spikes to be 2fc + 1 times the amplitude at the point of 
gradient catastrophe. 

2 A short review of the zero dispersion Umit of the inverse scat- 
tering transform 

Given an initial data (potential) q(x, 0, e) for the (1-1) that is decaying as x ±oo, the direct scattering 
transform for the NLS ([30]) produces the scattering data, namely: the reflection coefficient ro{z,e) and 
the points of discrete spectrum together with their norming constants (solitons). The time evolution 
of the scattering data is simple and well-known. Thus, to find the evolution of a given potential at a 
time t, one needs to solve the inverse scattering problem at this time t. Equivalently, we can stipulate 
that the initial data is assigned directly through the scattering data; thus, one can produce a solution 
to the NLS (1-1) by choosing some scattering data and solving the inverse scattering problem for i = 
(initial data) and for t > (evolution of the initial data) . The latter approach allows one to avoid solving 
the direct scattering problem and addressing many delicate issues associated with it (see [25] for more 
details). Since we are interested in studying the generic structure of solutions of the NLS near the point 
of gradient catastrophe, it will be convenient for us to define a solution to the NLS (1-1) by its scattering, 
not by its initial data. 

In considering the semiclassical limit of (1-1), one has to consider the semiclassical limit of the corre- 
sponding scattering transform. For the case of decaying potentials of type (1-2), this limit was discussed 
in [27] (inverse scattering) and [25] (direct scattering) , where it was shown to be a correspondence between 

a{x,t) = -^<i>^{x,t) + iA{x,t) (2-1) 

on the potential side and 

f{z) = foiz)-xz^2tz^ (2-2) 

on the scattering side. Here t > is fixed and fo{z) has the meaning of the "scaled" logarithm of the 
reflection coefficient ro(z,e) that corresponds to the initial data (1-2), that is, 

fo{z) ^ ^liraeroiz,e) . (2-3) 
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For example, the limiting fo{z) that corresponds to the potential -^^^^^^e ^ is given by (see 

[26]) 



/oW = (l-^) 



ln(l - z) 



+ zlnz + ln2 + -e, 



when 3z > 0. 



(2-4) 



Moreover, if x{a) denotes the inverse function to a(x,t), at t fixed, then x{a) and fa{z) are connected 
through a certain Abel type linear integral transform, see details in [25]. 

Since we have taken the perspective that we are starting with the scattering data of the form ro(z, e) := 
Q-irfoi^) we shall assume that ([27]): 

• fo{z) is analytic in the upper half-plane and has a continuous limit on M; 

• fo{z) ~ 0{z) as z — > oo; 

• there exists an interval (^_,/i+) C M such that 3/0(2:) > for z G (/i_,/i+) and 3/o(z) < for 

z e (-00, /i_) u (/i+,00); 

• 5/0 (z), z € ]R, has simple zeroes at ijl± and its values are separated from zero outside neighborhoods 
of /i±. 

Due to the Schwarz-symmetry of the scattering data for the NLS, /o(z) is also a Schwarz-symmetrical 
function, i.e., /o(z) — fo{z); generically, 3/o(z) has a jump across the real axis. The above assumptions 
are not overly restrictive: it was shown that in the case of a solitonless potential the formula (2-3) yields 
such an fo{z) after some process of analytic extension, or "folding" of certain logarithmic branch-cuts 
onto the real axis. In fact, for our goals it is sufficient to replace the analyticity requirement by a weaker 
condition: the contours of the RHP for the ^-function (which will be introduced in the next section) 
always lie within the domain of analyticity of fo{z) in C-^.. Generically, ro(z) has a jump discontinuity 
along the real axis due to the discontinuity in S/q. 

At any time t, the inverse scattering problem for (1-1) (in the solitonless case) with a fixed (not 
infinitesimal) e is reducible to the following matrix RHP. 

Problem 2.1 Find a matrix T{z) analytic in C \ M such that 



r+(z)-r_(z) 



roe 



rol^ + 1 roe 

1 



z e 



r(z) = 1 + -Pi + 0(2"^) , z-^oo. 



(2-5) 
(2-6) 



(In the case with solitons, there are additional jumps across small circles surrounding the points of discrete 
spectrum.) Then 

g(rz;,i,e) :=-2(Pi),2 (2-7) 



is the solution of the initial value problem (1-1) for the NLS equation. 
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The jump matrix for the RHP admits the factorization 



rp + 1 r ' 




■ 1 r ' 




■ 1 " 


r 1 




1 




r 1 



(2-8) 



where r = r{z;x,t) = ro(2;)e^(^*^^+'^^). 

Inspection of the RHP shows that the matrix r(z)* (where * stands for the complex-conjugated, 
transposed matrix) solves the same RHP with the jump matrix M(z) replaced by M~^(z) and hence 

Proposition 2.1 The solution r(z) of the RHP for NLS has the symmetry 

mimr ^ i . (2-9) 

In order to study the dispersionless limit e — > 0, the RHP (2-5)- (2-6) undergoes a sequence of transfor- 
mations (that are briefly recalled in Section 2.2) along the lines of the nonlinear steepest descent method 
[11, 26], which reduce it to an RHP that allows for an approximation by the so-called model RHP. The 
latter RHP has piece-wise constant jump matrices (parametrically dependent on a;,i,e) and, in general, 
can be solved explicitly in terms of the Riemann Theta functions, or, in simple cases, in terms of algebraic 
functions. The g-function, defined below, is the key element of such a reduction. 

2.1 The (yf-function 

Given /o(z), we introduce the g-function gi^z) = g(z;x,t) as the solution to the following scalar RHP: 

1. g{z) is analytic (in z) in C \ 7.,„ (including analyticity at oo); 

2. g{z) satisfies the jump condition 

g+ + g- = fo - xz - 2tz^ on 7^, (2-10) 

for a; e M and t > 0, and; 

3. g{z) has the endpoint behavior 

g{z) — 0{z — a) 2 -|- an analytic function in a vicinity of a. (2-11) 

Here: 

• 7m is a bounded Schwarz-symmetrical contour (called the main arc) with the endpoints a, a, ori- 
ented from a to a and intersecting M only at 

• g± denote the values of g on the positive (left) and negative (right) sides of 7^; 
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Figure 4: The typical zero-dispersion phase diagram for a one-hump initial data. Representative at 
different points in the (a;,t)-plane are the level-curves of ^{h). Only the upper-half spectral plane is 
depicted. The shape of the level curves was obtained numerically on a simple example. The plot of the 
amplitude of q{x,t,e) corresponds to the initial data q{x,0) = sech(a;), e = This is a pure-soliton 
case, which was used here only for the purpose of an effective illustration, as our results are valid for any 
generic (nondegenerate) gradient catastrophe. 

• the function /p = foiz), representing the initial scattering data, is Schwarz-symmetrical and Holder- 
continuous on 7„. 

Taking into the account Schwarz symmetry, it is clear that behavior of g(z) at both endpoints a and a 
should be the same. 

Assuming /o and 7™ are known, the solution g to the scalar RHP (2-10) without the endpoint condition 
(2-11) can be obtained by the Plemelj formula 




(2-12) 
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where R{z) = ^ {z — a){z — a). We fix the branch of i? by requiring that Xwnz^ao ^^^^ = 1. If /oC^) is 
analytic in some region S that contains 7™ \ tlie formula for g(z) can be rewritten as 

R{z) f /(C) 

where 7m C 5 is a negatively oriented loop around 7„i (which is "pinched" to 7™ in where / is not 
analytic) that does not contain z. Introducing function ft, = 2g — /, we obtain 

where z is inside the loop jm- The endpoint condition (2-11) can now be written as 

h{z) = 0{z — a)^ as z a, (2-15) 

or, equivalently, 

■^^"^^ -dC = . (2-16) 



(C - a)i?,(C)H 

The latter equation is known as a modulation equation. The function h plays a prominent role in this 
paper. Using the fact that the Cauchy operator for the RHP (2-10), (2-11) commutes with differentiation, 
we have 

where z is inside the loop 7m. 

In order to reduce the RHP (2-5)-(2-6) to the RHP with piece-wise jump matrices, called the model 
RHP, the signs of 3ft,(z) in the upper half-plane should satisfy the following conditions: 

• 3ft,(z) is negative on both sides of the contour (main arc) 7^; 

• there exists a continuous contour 7c (complementary arc) in C+ that connects a and so that 
'<sh{z) is positive along 7c. Since Qh{z) > on the interval (— oo,/x_), the point /i_ in 7c can be 
replaced by any other point of this interval, or by —00. 

Note that the first sign requirement, together with (2-10), imply that Qh{z) — along 7^. Since the 
signs of ^h{z) play an important role in the following discussion, we call by "sea" and "land" the regions 
in C-(-, where '<sh{z) is negative and positive respectively. In this language, the complementary arc 7c 
goes on "land" , whereas the main arc 7™ is a "bridge" or a "dam" , surrounded by the sea, see Fig. 5. 
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2.2 Reduction to the model RHP 

We start the transformation of the RHP (2-5)-(2-6) by deforming (preserving the orientation) the interval 
(—00, , which is a part of its jump contour, into some contour 7+ in the upper half-plane C+, such 
that /i+ G 7+. Let 7_ be the Schwarz symmetrical image of 7+. Using the factorization (2-8) the RHP 
(2-5)-(2-6) can be reduced to and equivalent one where: 

• the right factor of (2-8) is the jump matrix on 7+; 

• the left factor of (2-8) is the jump matrix on 7^; 

• the jump matrix on the remaining part of M is unchanged. 

It will be convenient for us to change the orientation of 7+ , which causes the change of sign in the off- 
diagonal entry of the corresponding jump matrix. On the interval {11^,00) we have 5/0(2;) < and it 
appears that the jump is exponentially close to the identity jump and hence it is possible to prove that 
it has no bearing on the leading order term of the solution (2-7) (as £ — 0: see [26] for the case when /o 
is a one-parameter family that contains (2-4) and [27] for the general case). Therefore the leading order 
contribution in (2-7) comes from the contour 7 = 7+07-. In the genus zero case, the contour 7 contains 
points a, a, which divide it into the main arc 7„i (contained between a and a, and the complementary 
arc 7c = 7 \ 7m- According to the sign requirements (2.1), the contour 7™ is uniquely determined as an 
arc of the level curve Qh{z) = (bridge) that connects /i+ and a, whereas 7c can be deformed arbitrarily 
"on the land". Because of the Schwarz symmetry 2.1, it is sufficient to consider 7 only in the upper 
half-plane, i.e., it is sufficient to consider 74.. 

Having found the branch-point a, the (7-function g{z) and the contour 7^, we introduce additional 
contours customarily called "lenses" that join a to fi+ on both sides of 7™ (and symmetrically down 
under). These lenses are to be chosen rather freely with the only condition that S/i must be negative 
along them (positive in C_). This condition is guaranteed by (2.1). 

The two spindle-shaped regions between j„i and the lenses are usually called upper/lower lips (relative 
to the orientation of 7^. At this point one introduces the auxiliary matrix-valued function Y{z) as follows 

g^9(z)o-3 outside the lips, 

1 -e-fH^-) 

1 

" 1 e-^h(z) 

1 



y(z) = e~^'»(°°''"^r(z) 



'+1 



(2-18) 



-3(2)0-3 



in the upper lip in 
in the lower lip in C-)-. 

The definition of Y{z) in C_ is done respecting the symmetry in Prop. 2.1, namely 

Y{z) = {Y{zy)-\ zeC_. (2-19) 
The jumps for the matrix Y(z) are reported in Fig. 5. 
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Figure 5: The jumps for the RHP for Y. The shaded region is where < (the "sea"). The blue 
contour is the main arc 7„j, the black contour in C+ is the complementary arc 7^ and the red contours 
are the lenses. The green circles show the boundaries of Aq, and A+. 
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The model RHP. In the limit e — > 0, according to the signs (2.1), the jump matrices on the comple- 
mentary arc 7c and on the lenses are approaching the identity matrix 1 exponentially fast. Removing 
these contours from the RHP for Y{z), we will have only one remaining contour 7^ with the constant 



jump matrix 



1 

-1 



on it. This is the model RHP. Calculating the (1,2) entry of the residue at 
infinity (see (2-7)) of the solution to the model RHP, one obtains the leading order term of the genus 
zero solution as follows ([26]) 

qo{x,t,e) = 5a(x,t)e?*(^'*'; (2-20) 

where 

^{x,t) = Ag{oo;x,t). (2-21) 

A direct calculation of g{oo;x,t) that uses (3-67) yields (1-3). 

To justify removing contours with exponentially small jump matrices, one has to calculate the error 
estimates coming from neighborhoods of points a, a and /i+ (fora and these neighborhoods are 
shown as green circles in Fig. 5). This is often done through the local parametrices. We shall consider 
the construction of the parametrices near the point as already done and known to the reader, see [26]. 
The only information that we need is that these parametrices allow to approximate the exact solution to 
within an error term £{z) — 1 -\- 0{e) uniformly on compact subsets of the genus zero region. 

3 Analysis near the gradient catastrophe point 

Let a = a{x,t) G C-|- be the branch-point in the genus zero region, where {x,t) is close to the point of 
gradient catastrophe {xo,to). 

For generic values of {x,t) the function h{z;x,t), according to (2-15), has the behavior ^h{z;x,t) ^ 
0{z~a)i; at the point {xo,to) of gradient catastrophe the behavior is instead ^h{z ; XqjIq) ^ 0{z — a)i. 
Thus for (x, t) in the vicinity of this point we obtain 

-h(z; x,t) = -(z- a)3/2 ^ cUz - a) + Oiz - af) , (3-1) 
e e 

where a = a{x,t) is the branch-point and Cq,Ci are some functions of x,t. The gradient (umbilic) 
catastrophe point is the one for which Co(a;ojio) = but Ci(0,0) 7^ 0, this latter inequality being our 
standing genericity assumption. 
According to (2-17) and (3-1), 

Co = Co{x,t) = hm — = — — / — dC . (3-2) 

z^a R{z) 3n (C - a)i?(C)+ 
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Lemma 3.1 The value of Ci at the point of gradient catastrophe is given by 

2V2fb f no 

Proof. To obtain Ci, we notice that at the point of gradient catastrophe {xo,to) 

where z is inside the loop 7„i. (This formula is not correct when (x,t) ^ (xc^o)-) Then, similarly to 
(3-2), 

Ci = Ci{x,t) = hm — = — — / — dC . 3-5) 

z^a R{z) IStt (C - a)i?(C) + 

Q.E.D. 

The goal of this section is that of introducing a suitable conformal coordinate C near z = a as in the 
definition below. 

Definition 3.1 (Scaling coordinate) The scaling coordinate (^(z) — (^(z;a;,t,£) and the explo- 
ration parameter r — T{x,t,£) are defined by 

-h{z;x,t)^U'^ +tC^, (3-6) 

£ O 

where C{a) = 0. 

The expression (3-6) is the normal form of the singularity defined by h{z; x, t) (in the sense of singularity 
theory [2]). 

The detailed analysis of r(a;,i;e) on space-time will be accomplished in Sect. 3.1; for the remainder 
of this section we dwell a bit on the details of the construction of r, C starting from the power-series 
expansion of h{z; x, t). 

Let us denote the expansion of h(z,x,t) as 

ih{z;x,t) = Co{x,t){z - a)'^ + Ci{x,t){z - a)^ + 0{z - a)i . (3-7) 
For At = Q = Ax we have 

-h{z;xo,to) = Ci(z-a)3 {l + 0{z-a)) (3-8) 
e 

and then the function C and the parameter t are defined by the formula 

2 

/ 5i \ ^ i 4 5 

Ciz) -.^ C{z;xo,to,e) := i—h{z;xo,to)] ^ ^h{z; Xo,to) = -(^ , r = 0. (3-9) 
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Thus, the function h{z; XQ,tQ) has a singularity (in the sense of singularity theory, i.e. the study of 
normal forms of degeneracies of critical values) at z = a. For x ^ xq, t ^ to this function undergoes 
a (smooth) deformation by which the coefficient Co{x,t) acquires a nonzero value that, consequently, is 
inherited by t. 

In the language of singularity theory this defines a (partial) unfolding of the singularity. It is a 
standard theorem [2] that for any such deformation there is a family of changes of coordinates z w so 
that 

ih{z] xo,to) = + T{x,t)w^ , w — w{z;x,t),w{a{x,t);x,t) = 0, (3-10) 
where T{x,t) and w{z;x,t) have the same smoothness class as the family of the deformation. 



Remark 3.1 To he more specific, changing variable from {z—a) to q = \J z — oi we then have a singularity 
for ih of type A4, with additional symmetry 

ih{q) = -ih{-q) . (3-11) 

Then the theorem guarantees the existence of a conformal change Q — Q{q) such that any deformation 
can be recast into 

ih{q; X, t)^Q^ + TiQ^ + TsQ^ + T^Q + To, (3-12) 

where Q is a local biholomorphic equivalence depending analytically on the deformations. The oddness 
forces To ~ T2 = and the fact that our particular deformation for ih starts with q^ forces Ti — 0. Since 
the theorem guarantees the existence of such analytic family of change of coordinates, a computation 
manipulating series allows to easily set up a recursive algorithmic procedure to find this function. See the 
next Remark 3.2. The most pertinent reference is Chapter 8 in [2]. 

Remark 3.2 It is not hard to find -recursively- the expansion of w and T in terms of the coefficients of 

the series of ih. If we set (for brevity we shift a to the origin, without loss of generality) 

00 00 
i/i(z) = CqZ^ + Cizs + ^ Cjz"T^ , w = wiz + ^ WjZ-' (3-13) 

we can equate the series expansions of 

ih{z) = wi + Twi . (3-14) 

From the coefficient of z^ we have T= and all the remaining coefficients of the expansion of w can 
be determined in terms of wi and the Cj 's. The first few are 

2 {w^^ - cj ) {wl + wl Cf + wiCf + Cl + Cl)wi ^ ^, 

6 Go 

{9wl - Sc^wl - 6C2C0 - Cf)w, 
ws = . (3-16) 
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The requirement that each term in the expansion should he analytic at Ca = Q determines Wi uniquely. 

2 

For example from (3-15) we must have Wi — +0{C{)): plugging this into (3-16) one sees that there 
can he at most a simple pole at Cq = and setting the residue to zero we determine the next coefficient 
in the expansion of wi . For example we have 

.Nc. + ^ + ^'^^^^--'^^'' +c^(q;). (3-17, 

25C^ 625Ci= 

^,^C!-^^+0{C^), (3-18) 

T = ^ + 0{C^). (3-19) 

While it is clear that this recursive procedure determines a formal expansion whose coefficients are analytic 
at Co = 0, it is not clear whether the expansion should he convergent. However the above-mentioned 
theorem guarantees the existence of such analytic expansion and hence it must coincide with this formal 
manipulation. 

In order to translate the normal form (3-10) into the desired one (3-6) we need to perform a simple 
rescaling 

w = el(^^\, T = ei(Jj%. (3-20) 

The function ^ is locally univalent in a neighborhood of z = a and C(q^) = 0- The function r is 
analytic in Cq at Cq = 0. Their local behaviour is 

^^\^ 6C0C2 , ,^,^2- 



C = eM-Cij i^l-^^+0{C^)]{z-a)il + 0{z-a)), (3-21) 



3 

4 



T = e-sCo(^— j (l + 0(Co)). (3-22) 

We point out that the determination of the root is fixed uniquely by the requirement that the image of 
the main arc (cut) where 3?/i = be mapped to the negative real ^-axis. 



Repeating identical considerations for the behaviour of h near a we define C(z; x, t) — (^(z; x, t). Before 
we can proceed with the detailed asymptotic analysis of the Riemann-Hilbert Problem 2.1, we need to 
establish more precisely the relation between the complex parameter r(a;,i,e) and the (x, t)-plane. 

We shall consider the scaling limit in which r is uniformly bounded; this means that T{x,t) in 
(3-10) must tend to zero as a 0(e5). This means that we will be considering some shrinking neighborhood 
of the point of umbilic catastrophe (2:0,^0)7 which will be determined in more details in Sect. 3.1. 

Remark 3.3 Incidentally, one could construct critical initial data for which there is a more degenerate 
gradient catastrophe -h{z; XQ,tQ) — -{z — a)^i^ {Ck+0{z — a)) with Ck 7^ 0. The case k = correspond 
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to a regular (non-gradient catastrophe) point (a;o,io); where the local parametrix is written in terms of 
Airy functions. The case k = \ is the one under scrutiny now and corresponds to a parametrix written 
in terms of Painleve I. For k > 2 it is easy to speculate that the PI parametrix needs to be substituted by 
a member of the Painleve I hierarchy. This will be investigated elsewhere. 

3.1 The map {x,t) h-)- T{x,t;e). 

The goal of this section is to determine the dependence of Cq,t, etc. on the space-time variables {x,t). 
Here and henceforth we use the notation ao — oq + ibo — a(xQ, to). 

Theorem 3.1 Near the point of gradient catastrophe (a;o,to) the behavior of a is 



Aa^ = K^{Ax + 2{ao + ao)At) + 0{At^ + Ax^ 



no 



-dc 



(3-23) 
(3-24) 



(C-ao)i?(C)- 

Proof. The branch-point a{x,t) is determined implicitly by the modulation equations (2-16) that can 
be written (see [26]) as 



F(a, a, X, t) — 



1 

2in 



f'{w)Sw f wf'(w)5w 

1^ R+iw) j^. 



= 0, 



R+{w) 

f{z; X, t) = /o(z) - zx~ 2tz'^, 



(3-25) 
(3-26) 



where 7m is a closed contour around the main arc and R{z) is chosen with the determination that behaves 
as z for 2 — oo. The Jacobian of F is ([27], Lemma 3.4) 



dF 



d{a, a) 



R{z) \z=a 
h\z) I 
^ R(z) \z=a 



h'(z 



R{z) \z=o 



a 



l(z) 
h'{z) \ 
R(z) I 



On account of the Schwartz symmetry h{z) — h{z), the determinant of this matrix is 



det 



dF _ ib 
d{a,a) ^ 2" 



h'{z) 



R{z) 



(3-27) 



(3-28) 



For {x,t) away from the gradient catastrophe in the genus zero region we have h'{z) — Co{z — q;)2 (1 + 
0{z — a) with Co — Co{x, t) ^ 0, so that the Jacobian (3-27) is invertible (as long as & = 5a > 0) and the 
standard implicit function theorem yields a{x^t) as a smooth function of {x^t). At the point of gradient 
catastrophe we have h{z) = Ci(a;o, io)(-2 — a)i(l + 0{z — a)) with Ci — Ci(xo,io) 7^ 0; therefore the 
matrix (3-27) is not invertible and -in fact- it is the zero matrix. 



Remark 3.4 Note that the Jacobian matrix -Jp-rs is 

d(x,t) 



dF 

d{x,t) 



1 

2iTr 



Sw 
Riw) 
w5w 
R{w) 



f R(w) 
T 4w Sw 
J Riw) 



1 4a 
a 4a2 - 252 



(3-29) 
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where a = a + ib. For any a € C+, we have 



det 



OF 

d{x,t) 



-26^ ^ . 



We expand each component Fj, j — 1,2, oi F around ao = ct{xa,to) and its complex conjuj 
denote only the dependence on a with the understanding that F depends also on a) as 



Fj {a; X, t) = Fj (ao; x, t) + d^Fj Aa + daFj + ^ [Aa, Aa] Hj 



Aa 
Aa 



+ 0(|Aan, 



where Aa = a — ao and Hj denotes the Hessian of Fj evaluated at (ag, x,t): 

w^f'{w)Sw f w^f'{w)Sw 



{w — a)^i?+(w) 
f'(w)6w 



(w — a)(w — a)R^{w) 
The off-diagonal entries of Hj vanish because 



{w — a){w ~ a)R-f-{w) 
f {w)Sw 



1 



1 



Stti J (w — a){w ~ a)i?+(w) 



IGnb 



f {w)6w 



(w — a)i?+(w) / {w — a)R+{w) 



f {w)6w 



1 



167r6 



,h'{z) 



R{z) 



= 



Let us denote by Gj^i the (1, 1) entry of Hj [j = 1, 2), so that 



(3-30) 
;ate (we 

(3-31) 



(3-32) 



(3-33) 
(3-34) 



3 



f {w)5w 
{w-aYR+{w) " 4 (5z V" ^(^) 



Z 5 f ,h'{z) 



J-0,1. 



Using the fact that g^r=y — 0, we then have (we suppress the (x, t) dependence for brevity) 

F,{a) = Fj(ao) + ^{G.^iAa') + 0(|Aa|3). 
Expanding also in x,t near a;o,toj we obtain 

= Fj{a) = Fj^^Ax + Fj^tAt + 3* (Gj^iAa^) + 0(1 Aap) + 0{At^ + Ax^), 



where Aa; — x — xq^ At = t — to. Equation (3-37) shows that Aa is of order y^|Axl~+~|At[. 
equations (3-37) (j — 1,2) for Aa using the expressions for the x, ^-derivatives in (3-29) yields 
2 



(3-35) 

(3-36) 

(3-37) 
Solving 



(Aa)^ 



[(Gi - aGo) Ax + (4aoGi - (4a^ - 2bl) Go) At] + 0{At^ + Ax^) (3-38) 



GqGi — GiGo 

At (xo,io,ao) we have Gi = agGo and hence a simplification of the above equation yields 



Aa" 



1 

Go 



[Aa; + 2 (ao + ao) At] + 0{At'^ + Ax^ 



(3-39) 



From (3-35) we find Gq = - with Gi defined by (3-3) and, thus, obtain (3-23). Q.E.D. 
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Corollary 3.1 Let us introduce the scaling variables x = xq + X, t = tg + e^T in the neighborhood D 
of{xo,ta). Then 



(X + 2K + ao)r)Ml + 0(£')), (3-41) 



(a-ao)2 =e5K2(X + 2(ao + ao)r)(l + C'(£5)), (3-40) 

_ Ai^fEb^ / 4 
~ 3K \JCi 

where K is given by (3-24) 0'''^d uq = ag + i6o- If^ the above formutE \/2iba stands for c^'^^^\/2ba (bo > 0), 
all the roots are principal and the argument of Ci is determined in such a way that the direction of the 
main arc is 

arg(z-ao) ^arg(Ci). (3-42) 
5 

so that the main arc is mapped to the negative real (-axis^. 
Moreover (3-40) can be written using (3-41, 3-24) 

Aa^ei^, (3-43) 

Proof. It is known from the Whitham modulation equations [27] that 

d^h{z; X, t) = R{z) , dth{z] x, t) = 2{z + a)R{z), (3-45) 



R{z) := — a)(z — a) , a — a{x, t) ~ a + ib . (3-46) 

Expanding in series near z — a and comparing the terms we find 

2iV2fb 4iV2fbia + g) 
CxCo = — , atCo = (3-47) 

Since ax = Aa^, using Thm. 3.1 we have 

„ 2iv2ib 2 4i\/2i6 , z — , 2,, 

Co = ^ = -£5 — ^ + 2{a + a)T (1 + 0(£5)) (3-48) 

6ax 3 /sV^ 



15iCi 

and the formula for t follows from the expression (3-22). Q.E.D. 

Corollary 3.1 states that -in the scaling limit- the map {xo + ei X, to+e^T) t'^ is a diffeomorphism. 
For definiteness and later purposes we introduce the following definition. 

Definition 3.2 The map 

v{x,t;e) = ^T\x,t;e) (3-49) 

will be called the Painleve coordinatization. 

^ _ 2 2 

°Since C, = e = C{z — ao){l -I- . . .) and C = (4/5Ci) = , the condition for the main arc is arg(Q = tt, whence the formula 
(3-42). 
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Figure 6: The level lines of the 5ft[i/i] in the C-plane and for different values of t. When cj) := arg(r) = ^ 
OT (f) = ^ the complementary arc or one of the rims of the lens is pinched and the solution enters the 
genus 2 region. Thus the genus region corresponds to ^ < arg(r) < 

According to the previous analysis, the function v{x,t;e) is a local diffeomorphism in the neighborhood 
D of the point of gradient catastrophe. To this end we formulate the following corollary. 

Corollary 3.2 In terms of the scaling coordinates X — ^7/? ; T = the function v reads, to the 

leading order as 

^ = l^^ = -'{^{^ + 2K + ao)T)(l + 0{e'i)). (3-50) 
Proof. A simple manipulation from the Def. 3.2 and eq. (3-41). Q.E.D. 

Remark 3.5 The complex-valued function v(x,t]e) is a diffeomorphism that maps the neighborhood D 
of XQ,tQ of size C(e5) onto a neighborhood V of the origin v — Q uniformly bounded (in e); in later 
sections v will play the role of independent variable for the Painleve I. 
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3.1.1 The image of the genus two region 

We can now find tlie opening of the sector S in the r^-plane (and hence v-plane as well) that is the image 
of the genus two part of the neighborhood D of the gradient catastrophe point. 

The critical value of -h{z;x,t) — ^0 + r^i is given by 2( + — 0, so the critical value (c of C is 



Cc = "|t. Thus 



i 4 s 3 

-h{zc;x,t) = -Q +tQ 
£ 5 



4 y ■ ■ V 4 

The breaking curves are determined implicitly by '^[ih{zcr', x, t)] — and hence by the condition 



v{t) = ^ (-3r)5 +5t(-3t)^ 



(-3)5 r5 



0. 



(3-51) 



(3-52) 



Care must be exercised due to the presence of the fractional powers: recall that the choice of conformal 
parameter C has been made so that the main arc is mapped to C < 0; the C image of the critical point, 
determines whether we are in the genus zero or two region as explained presently. 
The breaking curves correspond to the first directions where z^(t) — starting from the t > or 
-which is the same- 

= ^ ^ = arg{Q = - + —. (3-53) 



cos I - arg(Cc 



Thus the two arcs of the breaking curves correspond to two rays amongst the ones below 



arg(T) e 



r 27rfc 



,keZ 



(3-54) 



In order to explain which rays we need to choose we have to consider the topology of the level lines of 



for different values of r. Due to the scale invariance {C, i-> A^C, t i-> X^t, A > 



we can restrict ourselves to studying the argument of r only. We will consider t 



4 id 



The level lines 



of the real part of ih{z)/e in the (-plane for different values of (j) are plotted in Fig. 6. The transition 
between the genus and genus 2 regions occurs when the connectivity of the complementary arc and/or 
the rims of the lens needs to change. This happens for arg r = ^ when the complementary arc is pinched 
between the "sea" {diih < 0) or for arg(T) = when the main arc is about to break into two arcs. 

The above discussion about the directions of the breaking curves can be summarized in the following 
lemma. 

Lemma 3.2 The asymptotic image of the genus zero part of the region D around the point of gradient 
catastrophe {xo,to) under the diffeomorphism v = v{x,t,e) in the limit e — >■ is the sector 



arg(u) e 



An Yl-K 



(3-55) 



in the Painleve v-plane. 
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This is so because the argument of v is twice the argument of r and from the previous discussion. 

The complementary sector of aperture ^ is the asymptotic image of the genus-two region in the 
Painleve plane; in the following section we will describe the asymptotics of q(x,t,s) in terms of the 
tritronquee solution, which has -conjecturally- poles only in such a sector [12]. 

Remark 3.6 (The angle between the breaking curves in the (x, , i)— plane) Since we know now 
that the breaking curves correspond to the directions arg(u) = we can compute the angle at which 

the two breaking curves meet at the point of gradient catastrophe. Using (3-50), we calculate 

dxv — X, drv — 2>c{a + a), , k :— —i\j -p^ = —i | I ' (3-56) 

/t™ (C-ao)fl(C)+'^'^, 

Now we have v = x {X + -iaT + 2ibT) (1 + 0(e5)). 

The breaking curves correspond to arg(t;) — 27r/5,47r/5 in the Painleve plane. Thus, we have argLi = 
|7r — arg >f and aigL2 = |7r — sigx respectively. These values of aigL define the rays 

_ tan(argL,) 
^" 26-4atan(argL,)''' 

on the physical plane that are tangential to the breaking curves at the point of gradient catastrophe. So, 
the angle O in the {x,t)-plane between the breaking curves is 

tan(c^2) tan(tpi) 
= 2b-4atan(^.) - 2b-4atan(^,) 

^ ' -1 , tan(y2) tan(yi) ' ^ ' 

' 26— 4a tan(ijP2) 2b— 4a tan(v3i) 

where (pj — arg Lj , j — 1,2. 

Example 3.1 In the case when initial data (1-2) for the NLS (1-1) is A{x) — — secha; and S'{x) = 
— /ztanha;, where /i > 0, the point of gradient catastrophe was calculated to be {xq^Iq) — ^0, 2(i^-\-2) ) ' 
corresponding value a(xo,to) = «V7*+~2 and the slopes of the two breaking curves at the point of gradient 
catastrophe - 

±m = ±—^=L= ,6 3-59 

see [26]. Since the diffeomorphism v{x,t,e) asymptotically (as e ^ 0) maps breaking curves (of slope 
±m) onto the rays aigv = |7r ± |7r respectively, we have (all angle equations are mod 2tt) 

aigD^ = -^^==[l±2{a + a)m] = \±^tt, (3-60) 
Vl + m- 5 5 



®the expression for m is provided in Theorem 5.6, [26], however the expression for m in Theorem 1.1, [26], should be 
replaced by its inverse. 
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where the vector it = {l,m) and denote the derivative in the direction oflt. Thus, 

„, 2bm 2 _i 2bm 4 

tan h are >r = -TT, tan h are >r + tt = -tt. (3-61) 

1 + 4am 5 4am - 1 ^ 5 ^ ^ 

Equation (3-61) defines the slope of the breaking curves at the point of gradient catastrophe {xQ,to) in 
terms of a{xo,to) and Ci{xo,to). Let us show the slopes (3-61), found in [26], are consistent with Lemma 
3.2. Substitution of the slopes ±m from (3-59) into (3-61) yields 

TT TT 8 4 , ^ 

2 T g +a^g>^ = -7r± -TT. (3-62) 

We now use (3-3) to calculate x. Considering for simplicity the solitonless case fi>2, we obtain 

157r J^{C-a)R{C) + 



where 3/'(C) = §signC(l — X[-t,t]), T — \J \ 1, was calculated in [26], Sect. 6.4. (The choice of 

branch of R{z) in (3-3) and elsewhere in this paper is opposite to those used in [26]. That is why the 
sign of "j^j^ in (3-3) is opposite to the one that would have been calculated according to [26]). Direct 
calculation of the latter integral yields 



2iV2ib 



1 1 

{T-a)R(T) ^ {T + a)R{-T) 



(3-64) 



so after some algebra we get Ci = — ^ i, . Then argCi = f +27rfc for some /c G Z. Taking into account 

15(m+2)4 

(3-4-3), we obtain 

t: IT IT 2TTk ,„ 

arg-^-^+4-20^^- (^"^^^ 
Substitution of (3-65) into (3-62) shows that (3-62) holds with k = —1 (mod 5). Thus, slopes (3-61) from 
[26] are consistent with Lemma 3.2. We also conclude that 

Example 3.2 Using the same example as in Ex. 3.1 with pL ~ 2 (see Fig. 1) one can verify that the 
main arc has *^ = ^f§ = f^if ^'^'^ hence the nodes are in the directions arg(w„o(ie ~ Vp) — jq + kir. 
Looking at the orientation of the bisecant of the sector of poles of the tritronquee we see that -in this 
case- the nodes are aligned perpendicularly to the bisecant. 

3.2 The behavior of the phase near the point of gradient catastrophe 



26 



The genus zero (Whitham) approximation qo{x, t\ e) to the semiclassical 
solution q(x, t; e) is the leading approximation and it is valid uniformly 
in the "genus zero" region; its dependence on x, t is determined by the 
modulation (Whitham) equations [26]. 

These equations can actually be utilized to extend the definition 
of qo{x,t;e) beyond the genus zero region where -however- the actual 
solution q{x,t;e) will have a different behavior (typically of oscillatory 
nature), see, for example, [4], where go was extended beyond the break- 
ing curve. It will actually turn out that qo{x, t; s) can still be used in a 

neighborhood of the point of gradient catastrophe as a "reference" for Figure 7: A numerically accu- 
rate representation of the main 
describing the actual behavior of q{x,t;e). For this reason we briefly ^^^^ ^^^^^ ^^^^^ jj^^^g ^j- cj^^-j 

analyze qQ{x,t;e) near (xo,to)- Example 3.2 

In the genus zero region, the leading order approximation of the 

amplitude and the phase of q{x,t;e), according to (2-20), (2-21), are given by b{x,t) and 4(jr(oo;x,i) 

respectively, where the branchpoint a{x,t) — a{x,t) -f ib{x,t) and the g-function g was defined by the 

scalar RHP (2-10), (2-11). In ([26] Lemma 4.3, formula (4.43)) it was shown that 




gj;{z;x,t) ^ ^ (^^/{z - a){z - a) ~ , gt{z;x,t) = (z + a)^y{z - a){z - a) ~ z^ 



(3-67) 



where the determination of the square root is such that they behave like z at infinity^. Hence for the 
phase <l>(x, t) — 4:g{oo; x, t) we have 



$^ = -2a{x,t) = -23?(a), 
$t = -Aa{x,tf +2b{x,tf 



-25R(a(a-f a)). 



(3-68) 



Theorem 3.2 The phase A$(x,t) :— <l>(x,<) — <I'(xo,to) near the point of gradient catastrophe (xo,to) 
has the expansion 



A$(x,t) =-23fJ|^ao [Ax+{ao + aa)At] j - 
= - 2aoAx - 2(2ag - bl)At - es^R 



— (Ax + 2(ao + ao)At)5 



2i t3 
CbY 



0{Ax^ + At^), 



0{Ax^ + Ai^) ^-69) 



(3-70) 



where Ax — x — Xq, At — t — to and K, C and r were defined in (3-24), (3-43) and in (3-22, 3-41) 
respectively. 

Proof. If wc write oi — + Aa we have from (3-68) 



= -2ao - 2Aa 



$t = -Aal + 2bl - 8aoAa + 46oA6 - 4Aa^ - 2A6^ 



(3-71) 



''Note that in [26] the determination being used is the opposite one. 
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In our problem Aa = 0{e^) and hence we can approximate 

$^ = -2^ao - 23?Aa , = -2^{ao{ao + ao)) - 23? (2(ao + oq) Aa) + 0{s^] 

From (3-40) and a straightforward integration we obtain (3-69).Q.E.D. 



(3-72) 



Remark 3.7 The formul<£ for $ allow us to extend the definition of qo{x,t) within the genus-2 region 
using (2-20); taking the imaginary part of (3-4-3) we find that 



h{x,t)=h^ + \ei^[y){l + 0{e^) 



and hence 

qo{x,t) = [ uo 



1 Ic-/"^ 



(3-73) 



(3-74) 



X exp - 
e 



$0 - 2aoA2: - 2(2ao - bl)At - e^'Si 




(1 + 0{Ax^ + At^)) (3-75) 



In the following we will understand that q(){x,t), b{x,t), a{x,t) have been extended as indicated above. It 
is to be noticed that this extension is discontinuous due to the definition of t (3-4-1) involving a square 
root. Such ambiguity will not be present in the final formul(B. 

4 The Riemann— Hilbert problem for Painleve I 

The heart of the present paper is in the detailed 
analysis of the "local parametrix". This will be 
constructed in terms of the so-called Psi-function 
^(^, v) of the Painleve I Lax system, that depends 
on the spectral variable ^ and the Painleve vari- 
able V. The analysis of the Ricmann-Hilbert prob- 
lem for 4'(^, v) is contained in a number of papers 
and books, sec, for example, [20, 13]; this analysis, 
however, does not cover the case when the Painleve 
variable v is at or is approaching a pole v = Vp of 
the solution to PI that is defined through 4'(^, t;) 
(as the isomonodromy condition). Furthermore, 
it can be shown ([22]) that '^{£,,v) has a pole at 
V — Vp. Analysis of the RHP for 'I'(^,w) at or 
close to a pole v — Vp of the tritronquee solution 
(transcendent) y{v) to PI is a matter of crucial 




Figure 8: The jump matrices for the Painleve 1 RHP: 



here 
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importance in our study of the height and the shape 
of the spikes. We start from the summary of the 

known facts about PI. Let the invertible matrix-function P = P(C,w) be analytic in each sector of 
the complex ^-plane shown on Fig. 8 and satisfy the multiplicative jump conditions along the oriented 
boundary of each sector with jump matrices shown on Fig. 8. The entries of the jump matrices satisfy 
the following symmetry conditions 

l + /3o/3-i = -^2, (4-1) 
l + /3_2/3_i =/3i, 

so that the jump matrices in Fig. 8 depend, in fact, only on 2 complex parameters (that uniquely define 
a solution to PI). The matrix function P(C,w) is uniquely defined by the following RHP. 

Problem 4.1 (Painleve 1 RHP [20]) The matrix P(^;w) is locally bounded, admits boundary values 
on the rays shown in Fig. 8 and satisfies 



P+ = 



P_Af, 

^CT3/4 



/ + o(r') , 



(4-2) 
(4-3) 



where the jump matrices M — M(^; v) are the matrices indicated on the corresponding ray in Fig. 8. 

For any fixed values of the parameters /S^, Problem 4.1 admits a unique solution for generic values of 
V] there are isolated points in the u-plane where the solvability of the problem fails as stated and it will 
need to be modified. 

The piecewise analytic function 5'(^,i;) — P(^, i;)e''°"=', where "i? :— t?(C;w) — — vS,^ , solves a 
slightly different RHP with constant jumps on the same rays. The new jump matrices can be obtained 
from the old ones by replacing the exponential factor in every jump matrix by one. It then follows that 
it solves the ODE [20] 

d 



y' 2e + 2y£,-v + 2y2 
2e - 2y -y' 



(4-4) 



where y = y(v) solves the Painleve I equation 



y 



6y' 



(4-5) 



Direct computations using the ODE 4-4 and formal algebraic manipulations of series along the lines 
of [29, 17, 15, 16] show that ^' admits the following formal solution 



V2 



1 -i 
1 i 

~Hia3 ~Hjl + ya2 , iv^-iH]-2y') iy'-2tHiy 
+ 24^1 + 4Ci 



J173 
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H, 



1 



yv - 2y3 



(4-6) 



where 0(^^^) denotes the sum of terms with higher order powers of Such an expansion has to 

be understood as representing the asymptotic behavior of an actual solution of the ODE (4-4) within a 
sector of angular width smaller than ^ . 

4.1 Failure of the Problem 4.1 

The choice of the parameters is (transcendentally) equivalent to the choice of Cauchy-initial values 
for the ODE (4-5) ; it is known since the original work of Painleve that the only (finite) singularities of 
Eq. (4-5) are poles and these poles coincide precisely with the set of exceptional values of v for which 
Problem 4.1 fails to admit a solution. From the PI equation (4-5) for y{v) one can find the Laurent 
expansion around any such pole v = Vp to be of the form 

yiv) 



(4-7) 



We can then proceed as follows [22]: define the matrix ^'(^,i;) via 



$(^;f) 



-T3/2 



2(?-!/) 

1 





1 -K^' 



(4-8) 



1 



2(?-!/) 



T3/2 



It then satisfies the ODE 



Vi^;v) ■.= 2e 



~ 2y^ + yv + -{y 



l\2 



2 
V{i-v) 

y' 



$(^;^;)e''^^ 
3 



2(e-y) m-yf 

It is promptly seen from a direct computation that the function 1^(^; v) admits a limit as v 

=^H^ 

' ^ ^ y'i , 3 



(4-9) 
(4-10) 



= 2f -v(, + Hi 



- 2r + vy 

y'£, 



y{i-y) m-y? 



y_ 

2y^yi(-y)^ 8{^-yy ' 
^ 2e - - 14/3, 



Hr 



H, 



(4-11) 
(4-12) 



y_ 

2y 6 

where the convergence is uniform over compact subsets of the f plane [22] . It was also shown ibidem that 
^'(^;'y) tends to a finite (holomorphic) matrix ^'(^;Wp) which satisfies the (essentially a scalar ODE) 



2 

- - 14/3 



^'(^;wp) =:^(^;i;p,/3)*(C;Wp). 



(4-13) 
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Most importantly, the solutions to the system (4-9) and ^(^;Wp) to the limiting system (4-13) 

have the same Stokes' matrices. In fact, the Stokes' matrices for these solutions are the same as those 
of Vl'(^, v), except minor changes introduced by the obvious nontrivial monodromy of the transformation 
G{£,,v) in (4-8). That follows from the isomonodromic property of the equation (4-4) that defines the 
PI equation and the fact that the left multiplication by G(^, v) does not change the Stokes' phenomenon. 

Remark 4.1 The formal monodromy around ^ = oo for 4* is —ia2 but the one o/\E' is ia2 because of the 
additional monodromy (—1) around ^ = y. Using the explicit expression (4-8) the reader can also verify 
that 

(/ + 0(r^))e^"^ (4-14) 



1 i 
1 -i 



V~2 

4.2 Analysis in a neighborhood of the pole of PI 

It is essential for our application to investigate the behavior in which t; — > at a certain rate, namely, 
to study how (and in which sense) the limiting expansion (4-16) is approached. It is proven in Appendix 
A that 



^(?,i^) = r^"^'^(^i + 'T3) 



1 

i 



or 




where the term 0{^~^^^) is uniform w.r.t. w in a finite neighborhood of the pole v — Vp. 
The limiting case of (4-15) ior y — oo (i.e. for v = Vp) is 



(4-15) 



*(?;fp) 



V2 



-i -1 
-i 1 



14/3 



, 98/321 



4802 




(4-16) 



where the expressions for the various coefficients are obtained from the formal solution of the ODE (4-13) 
using the standard techniques in [29]. 

4.2.1 The tritronquee transcendent 

The term tritronquee dates back to Boutroux [6, 7]. A generic solution to the ODE (4-5) has infinitely 
many poles that accumulate asymptotically for large \v\ along the rays arg(w) = Certain one- 

parameter families (corresponding to the vanishing of one of the Stokes' parameters jSk of the associated 
Riemann-Hilbert problem (4.1) have the properties that along one of these rays the poles eventually stop 
appearing as — ^ oo, or they get truncated, whence the term tronquee. If two consecutive /3j's vanish 
we have 5 very special solution for which the poles truncate along three consecutive rays, whence the 
naming tritronquee. In fact there are -strictly speaking- several tritronquee solutions: they correspond 
to the vanishing of the /3j's on two consecutive rays in Fig. 8. There are -thus- 5 such functions. 
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However a closer look [20] reveals that the solutions y{v;{l3j}) have 
the symmetry 




yiv; {Pk}) = e^y (e'-^v; {/?fc+2„}) , n e Z , (5^+5 := /3fc (4-17) 



and hence there is essentially only one tritronquee solution. 



Such a solution is characterized by the following theorem. 



Theorem 4.1 ([20], Thm. 2.1 and Corollary 2.5 and eqs. (2.72)) 



There exists a unique solution y{v) corresponding to = Q — 



with the asymptotics 




Figure 9: The plane of the 
tritronqee transcendent, with 
the poles covered by disks of 
radius e^; given that the v- 
plane is an image of a es-size 
neighborhood D of the gradient 




It is conjectured in [12] that the tritronquee solution has actually catastrophe pomt, the transver- 



Whatever the standing of the conjecture, it is quite essential for us that the poles "want to be" in 
the sector of the v plane that indeed corresponds to the genus-two region under the identification of the 
z)-plane and the (scaled) X,T plane given by Def. 3.2. 

We are going to see below that -in fact- each pole of the tritronquee corresponds to a "spike" in 
the asymptotic solution of NLS and such spikes are to be expected only in the region of paroxysmal 
oscillations (genus 2). 

While our analysis does not rely in the least on the position of such poles, the "physical intuition" 
strongly suggests that indeed they will be confined to the indicated wedge. 

5 Leading order approximation of NLS away from a spike 

As we have seen in Section 3.1 (Def. 3.2), the map v = v(x,t;e) maps diffeomorphically a neighborhood 
D of size 0{e^) to the complex-^ plane, namely, to the plane of the independent variable of the Painleve 
I tritronquee transcendent. 

*Note that in [20] the independent variable x coincides with our —v. 




sal size of the spikes in the am- 
plitude oscillations of q{x,t;e] 
is 0(e). 
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As we have seen in the previous section, the region of the {x, t)-plane that corresponds to the genus-0 
region is mapped to the complement of the sector arg(w) e (27r/5, 47r/5). In the present section and the 
following Section 6 we shall consider two different limits in which a point {x,t) approaches the point of 
gradient catastrophe (a;o,io)- These limits, express in terms of the diffeomorphism v{x,t;e), are: 

• v{x,t;e) is in a compact subset of the "swiss-cheese" region Ks, S > is constant^, so that v is at 
least on the distance S away from any pole of the tritronquee solution; this is the case considered 
in this section; 

• v{x,t;e) G Bs, where Bs is a disk of radius S — 0{ei^'^), v > 0, centered a,t v — Vp - one of the 
poles of the tritronquee i.e., v can approach a pole Vp of the tritronquee at a rate es or faster; this 
is the case considered in Section 6 . 

Of course, we shall consider also the case where v is exactly at a pole (Section 6), as well as the case 
when v{x,t;£) — ^ Vp at the rate 0{ei~'^), where v G (0, ^) (this section). What will transpire from the 
analysis is the following enticing picture: 

• The spikes are in one-to-one correspondence with the poles of the tritronquee solution; 

• The transversal size of the spikes is of order 0{e); 

• The separation of the spikes is of order ©(et); 

• The amplitude of the spikes is three times the amplitude of qo{xo,to), i.e., the amplitude 
predicted by the modulation equations for the genus-zero region; 

• -Each spike has the shape of a rational breather. 

5.1 Asymptotic behaviour away from the spikes 

In the genus zero region, the leading order solution to the RHP for Y{z), i.e., solution to the model RHP, 
is ([261) 



1 


— i 


-1 ■ 




i 


1—1 


2 


1 


i 


\z-a) 


-1 





*o(2) = 



More than its specific form, it is important that near the point z ~ a it has the behavior 



(5-1) 



1 



vI'o(C)=0(z-a)-^(z-«)T 



i 1 
-1 -i 



(5-2) 



with the jump matrix i(T2 



1 

-1 



on the main arc. Here ©(C — o) denotes a matrix function that 



is invertible and analytic in a neighborhood of a. 



'The definition of is in Sect. 1.1. 
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We shall construct an approximation to the matrix Y{z; x,t, s) appearing in (2-18) (and henceforth 
to the matrix F) in the form 

{£{z)'i>o{z) for z outside of the disks DcD^ 

£iz)'^o{z)Va{z) for z inside of the disk D„, (5-3) 
£{z)'iQ{z)Vaiz) for z inside of the disk D^, 

where Aq, are small disks centered in a, a respectively, see Fig. 5. 

Remark 5.1 There is an additional parametrix for Y in a disk A-^- centered at the point z = ; however, 
as shown in [26], this parametrix has a uniform estimate 0{e), so it will not affect the accuracy of any 
of our expressions. 

Due to the symmetry of the problem in Prop. 2.1 we must have 

r^{z) = {r^{z)r-' (5-4) 

and hence it suffices to consider the construction near the point a only. 
5.1.1 Local parametrix 

The local parametrix Va{z) = Viz) (we understand and suppress the subscript a) must satisfy a certain 
number of properties (see Thm. 5.1), one of them being the restriction of V{z) 



Viz) 



= 1 + 0,(1) (5-5) 



on the boundary of Dq,, where Oj(l) denotes some infinitesimal of e, uniformly in z G dDa and in 
{x,t) £ K ~ v^^{Ks) (here a small (5 > is fixed). 

If V (and the corresponding parametrix near z = a) can be found that satisfy those requirements 
then the "error matrix" £(z) is seen to satisfy a small-norms RHP and be uniformly close to the identity. 
More precisely the matrix £ has jumps on: (a) the parts of the lenses and of the complementary arcs that 
lie outside of the disks Da, Da, and; (b) on the boundaries of the two disks Dq,Dq. The jumps in (a) 
are exponentially close to the identity jump in any norm (including L°°) while on the boundary of 
the disk D^, we have 



£+{z)=£^{z)^o{z)V{z)%\z) 



= (1 + 0,(1)) (5-6) 



Since £{z) = l + 0{z~^) as z — oo it follows [11] that ||£(z) — 1|| (uniformly on the Riemann-sphere) 
and that the rate of convergence is estimated as the same as the 0,(1) that appears in (5-5) as e 0. 

Thus, the accuracy of the approximation (i.e. neglecting the term £{z)) is directly related to the rate 
of convergence to the identity matrix of the local parametrix V on the boundary of the disk(s). 
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Definition 5.1 (Local parametrix away from the spikes.) Let C{z;e) be the local conformal coor- 
dinate introduced in Def. 3.1 so that 



-h{z;x,t)^e{C;T) = U~^+tC-- 
£ 5 



(5-7) 



Let ^'(^;t)) denote the Psi-function of the Painleve I problem with /3q = = (and /?_2 = ^1,. 

/32 = — Ij- The parametrix V{z) is defined by 



Viz) 



1 

7^ 



1 1 



1 

-1 



,0(C;t)£T3 



(5-8) 



Theorem 5.1 The matrix V satisfies: 



1. Within fi> a, the matrix V{z) solves the exact jump conditions on the lenses and on the complementary 
arc; 



2. On the main arc (cut) 'P{z) satisfies 



V+{z) = a2V-{z)a2 , 



(5-9) 



so that '^o'P within solves the exact jumps on all arcs contained therein (the left-multiplier in 
the jump (5-9) cancels against the jump of^^^); 

3. The product "^q{z)V{z) (and its inverse) are -as functions of z- bounded within Dq,, namely the 
matrix V{z) cancels the growth of at z — a; 



4-. The restriction ofV{z) on the boundary of Da is 



V{z) 
where 



= 1 



Hj 



CT3 



16; ^ 2c 



H, 



16 



0-2 



+ 0(0, 



V — -T^, 

8 ' 



Hi = ]^{v'f +yv^2y^ ^ J y{s)5s. 



(5-10) 



(5-11) 



Proof. (1) The matrix '^{^; v) has constant jumps of the same triangularity as the jumps indicated 
in Fig. 8 (with /3o = = and /3_2 = = ^ 1 = —f3i)- These are altogether of the opposite 
triangularities we need, hence the second-last (constant) matrix in (5-8). The last multiplication with 
ge(C;t) _ gj/i(z;x,t) gjygg ^jjg exact (non- const ant) jumps on the parts of the complementary /main arcs 
and lenses within the disk D„. On the other hand, the matrix 



J-(z) := 



1 



1 1 



(5-12) 



has the jump T-i^{z) = ia2J^-{z) on the left, whence the part (2). 
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As for part (3), the product (2)7^(2) is a bounded function of z because the singularities of 5*0 (2) 
are canceled by those of 



1 


' -i -1 ■ 




' i 1 ' 


1 


— i z 


2 


1 i 


\z ~ a ) 


-1 -i 


7i 


1 1 



—I I 

1 1 



z — a 
z — a 



0{1) 



(5-13) 



since C(z) = 0(z — a). In fact, we see that the product is actually analytic. 

Finally, part (4) follows from the asymptotics of ^{£,; v). Indeed, for z e dDa the conformal coordinate 
(■(2) grows (homothetically) as e^^ and hence we can use the expansion (4-6) for ^(^;i;) near infinity. 
To see how it works let us recall the notation 

4.5 , 3 



(5-14) 



so that we can write 
Viz) 



1 


—i i 


/■ "3 T\°-3/4 


■ 1 -i' 


2 


1 1 


^ ^ ('^+2) 


1 i 



1 P=== + + 0{C^) X 



2C + T 



1 
-1 



=(5-15) 



1 1 



c 



0-3/4 



i 1 
-i 1 



yC+I 2C + r 



8C 



+ o{c 



vC 2C y 



We also have 

^(C;r)-^(^C+^;-^r^ 

so that -continuing from (5-16)- we have 



4,5 ,3 4 



+ 0(0, (5-17) 



7^(2) 



Hr 



(T3 1 

16; 7^ + 2c 



= 1 + Uf/ + 



\ 0-3 1 



16; vc 2c 



^6 



256 8 
16 



0-2 



+ 0(C-5) =(5-18) 



2/ 



0-2 



+ o(r^). 



(5-19) 



Q.E.D. 

At this point we already know that the error term £{z) in (5-3) is within 0{e^) from the identity; 
if we simply ignore it, we would get the leading order approximation to Y and -in turn- to F, which 
would produce the leading order approximation of the NFS solution q{x, t, e). Next, we will find the first 
subleading approximation by solving the first step in the iterative approximation of the error term itself. 
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5.2 Subleading correction 



Theorem 5.2 The behavior of the focusing NLS in the domain D of the point of gradient catastrophe 
(scaled like ei ) is given by 



qix,t,e) = 60 



1 - 2£53 



Cbo 



X exp 



$0-2 {oqAx + {2al - bl)At) + 2e^^ 



2i 

Cbo 



Hi{v) 



(5-20) 



where ©(e^) term is uniform "away from spikes", i.e., is uniform in (x,t) G Ks = ^ ^{^s) with some 



fixed S > 0. Here $0 = ^{^Q, to), ct{xo, to) = oq + ibo, C 

1 



'5C 



-) ' given by (3-44), 



Hi^^-{v'{v)f + vy{v)-2y\v) 



(5-21) 



is the Hamiltonian of the Painleve I equation, evaluated along the tritronquee solution appearing in 
Theorem. 4-1, while v can be expressed (Corollary 3.2) as 



-le 



2ibo 
1^ 



Ax -f 2(ao + ao)At (1 + ©(e^)) 



(5-22) 



uniformly in D. 



Remark 5.2 (Comparison with the Conjecture from [12]) The approximation formula (5-20) "away 
from the spikes" is consistent with the conjecture from [12] about the behavior of the amplitude and 
the phase of the genus zero (modulated plane wave) approximation qo{x,t,e) in the genus zero (non- 
oscillatory) part of the neighborhood D of the point of gradient catastrophe. This conjecture can be 
written as 

U + i^/lToV = Uo+ ^VC^Vb + e^Ky(v) + ©(et ), (5-23) 

where U — jgp, V = ed^ a.rg(q), with K being a (complex) constant. Using our expression (5-20), 
Uq = bg, and the fact that Hj{v) = y{v) we find 



boyjv) 
C 



+ ©(es), 



V = -2ao + 4£5 3fi 



yjv) 

C 



so that 



U + iboV = bl-Ae^bo^ 



c 



— 2iaobo + 4ie5 6o3? 



yjv) 

C 



Aibo 



bi - 2iaabo + £5 -^y{v) + ©(e^ ). 



(5-24) 

(5-25) 
(5-26) 



To replace the error estimate ©(£=) in (5-25) with the estimate 0{e'^) from the conjecture (5-23), if 
correct, would require calculation of the higher order corrections to the RHP (2-18). It may be true that 
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the approximation of q is in powers of rather than in powers of , but that, again, would require 
additional analysis. The situation here may resemble the analogous statements in random matrix theory 
in regards to the expansion of the partition function in even powers 1/N . We also did not dwell into the 
notation of [12] to compare all the constants used. Finally, we omitted the term proportional to (t — tg)^ 
in (5-23) because, in our scaling, it is of order and hence not "visible" at this order of approximation. 



Proof of Thm. 5.2. In a RHP of the form 

E{z) = l + 0(z-i) 



£+{z)=£_{z)[l + ^M{z) 
the solution (if it exists) can be written as 



£_(s)AM(s)ds 



(5-27) 
(5-28) 

(5-29) 



where the integration extends over all contours supporting the jumps (it is a simple exercise using 
Sokhotskii-Plemelji formula to verify that the above singular integral equation is equivalent to the 
Riemann-Hilbert formulation). If -in addition- the term AM(z) is sufficiently small in the appro- 
priate norms (at least in L°° and of the contours) then the above formula can be used in an iterative 
approximation approach, where 



1 

2iTT 



£'^\s)AM{s)ds 



J =0,1, 



(5-30) 
(5-31) 



s — z 

which can be shown to converge to the desired solution. 

In the RHP for £{z), stated in Section 5.1.1, we have AM exponentially small (in e) in any norm 
on all parts of the contour outside the disks Aq and Ag, and approaching zero in the L°° norm on the 
boundaries dDa, dDa. The latter estimate is valid in any norm due to compactness. We shall thus 
find the first correction term in the above approximation procedure. 

As noted in Theorem 5.1, part 4, and in (5-6), the jump of £ on SD^ is (note that commutes with 

0-2) 



£+iz) ^£^{z)^oiz)V-Hz)%\z), 



where 





t3 




- Hr + — 






16 




1 












"2C 







1 i 


1 


■ 1 -i ' 




i -1 




-i -1 


) 



0-2 



(5-32) 

(5-33) 
(5-34) 
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2 

P ■= - 



ui := Hi 



C VpC 



Ni := 



1 

1 -i 
-i 1 



u^l - U2a2) 



16 



Due to the symmetry, the jump on the disk around z = a is 




2Ct 



C^(z) :=C(^) , p"^ =p{z) 



(5-35) 
(5-36) 
(5-37) 

(5-38) 
(5-39) 



We now proceed to the computation of the first-order correction to 8{z) according to the formula 
(5-31). In that formula the integral should extend to all the jumps of £, which include the lenses, the 
complementary arcs and the disk around /i^ ; the former are exponentially small and the latter is of order 
0(e), thus they can be neglected altogether to within this order. 

When z is outside the disks this residue computation annihilates the analytic term with \Jp/Q in 
(5-35) and yields 



£(i)(z) 

( u\l- U2<72 



7il 



1 N J) t ~ a dt Ui f t — a 



dt 



(5-40) 



1 



St 



Ui + U2<J2 



\z-c.\=sC{t){t-z) 2m 



1 



6t 



+ 0(£5) =(5-41) 



_^ ^ e^^uxNx I a — a eiuiN^ I a ^ a 2 / u^l — U2<T2 ^ ufl + U2(J2 



2(z-a)V C 2{z-a) 



C 



\2C{z-a) 2C{z~a) 



(5-42) 



So, 



^""'^ ~2\~C z-a ~ Y\l z-a~^' \ 2C{z~a) 2C{z - a) 



+ 0(£5),(5-43) 



where -by definition- C appears as C(z) ~ e^C{z — a){l + ...). We now need to use once more (5-31) 
with the expression (5-43 and retain only the terms up to ei): 



^ ^ '-^ ^ 4 \C{z^a) ' C{z-a) 
Hence the approximation of the solution Y{z) is 

Y{z){z)^£{z)Mz)=n^)l 



ei ( u\NiNt ujN^NA i\ui\^ei fN^Ni NiN^ 



A\C\ 



z — a z — a 



—i 


-1 ■ 




i 


1 


1 


i 


\z — aj 


-1 


— i 



(5-44) 



(5-45) 
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Writing £{z) = 1 + ^ + 0{z ^) near z ~ od, we have 



Y = £{z) 1 



4z 



-i 

1 



-b 
b 



-Oiz-'). 



z 2z 

The correction comes from the (1, 2) entry of £i; we use the second iteration S^^^^ and 
Mies hib uies /-2i& 2 /w^l - 7/20-2 ull + U2a-2 



c 



4\C\ 



I 

V =4ct2 



1 Ml 2ib 1 Ml —2ib 

-ie ^ — \ ie^ — -i / 

2 V C 2 V C 



-"2 



/ "-Z \ 2 ZM? 2 iU? 2 |mi 



12 
2 



2C 2C 



\C\ 



. ,2ib \ £5 /2iu? 2iMi 4|mi|2\ 2„ /M2\ 

«e53fj ( J— Ml I + — ( ^ - + ^7^;^ ) +£53 l^-f j 



4 V C- 



We thus have 

q{x,t,£) = -2 Hm 2ri,2(;j) = _2eT'»(°°'^'*) hm z{£^n)„ = 
I 6(x,t) + 2z6£5 3fi ( 1 -2&£S5R 



\C\ 
~2i 



(^) 



7UI 



- 2£55 



C 



: exp 



+ 2i£53fJ I W — mi 




5(x,i)-2£i3('^ ) +0(£S) 



Now care must be exercised before factoring b(x,t) out: indeed from (3-43) we see that 
and thus the r in (5-51) cancels out and we obtain 



q{x, t, e) = bo exp 
= 60 exp 
= bo exp 



-<i>{x,t) + 2ie^^ 



:<^{x,t) + 2ie^^ 



'-^{x,t) +i£53fj 



Cb 
'~2i 



Ul 



l-2sl^l^+Oiel 



cb\"' + h) 



l-2£tc,(A)+0(£t) 



(5-46) 
(5-47) 



(5-48) 
(5-49) 

(5-50) 



2i 

CbY 



(5-51) 

(5-52) 

(5-53) 
(5-54) 
-0(£5)(l5.-55) 
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Note that we could also replace the remaining occurrences of b by bo in (5-55) since it would affect the 
result by terms of order ©(et). Recall here the expression (3-69) for the increment of phase near the 
point of gradient catastrophe: 

A$ = -2aoA2; - 2(20^ - bl)At - e^-'Si I \ + 0{Ax^ + Af ) (5-56) 

\ V Co 8 / 

The reader may notice that the discontinuous term containing cancels and we complete the proof. 
Q.E.D. 

Remark 5.3 From the expression (5-20) it is clear that the approximation cannot hold in proximity of 
a pole of the tritronquee y{v), for the Hamiltonian Hi has a simple pole there and y has a double pole. 

The reader could verify that the above analysis holds as long as we approach a pole v = Vp but not too 
quickly; 

v~vp^ Oiei-") , \>iy>0 (5-57) 

In this case the formula in the above theorem is still correct but with the error term of order 0{e^^^'^), 
and then the leading correction has -in fact- order 0(e^^'^). Indeed the term 0{C,^^) contains terms 
with triple poles (see (4-6)), and -in general- the term has a coefficient with a pole of order k in the 
Painleve variable v. Hence the estimate 0{e^ in (5-34 and following would have to be replaced throughout 
by 0(et-3>^). 

It appears that something awry is occurring when we approach a pole too fast, and a different 
approximation parametrix need to be constructed. This is the goal for the rest of the paper. 

6 Approximation near a spike/pole of the tritronquee 

With the preparatory material covered in Section 4 we shall now address the approximation of q{x,t;e) 
near a spike or -which is essentially the same- in a (shrinking) neighborhood of a pole of the tritronquee 
solution. 

In order to motivate the construction used below we illustrate the difficulties in constructing the 
leading approximation to the solution Y{z) of the RHP: it should appear that the local parametrix V 
needs to be expressed now in terms of the modified Psi function "if (4-8). 

Looking at the asymptotic expansion for large ^ of ^ (4-15) it appears that the first modification we 
need to make in order to match the boundary behaviour of the ^E* with the outer parametrix is to replace 
the solution '^q to the model RHP by 



^i{z) := 



1 


—i 


-1 ■ 




i 


1 


2 


1 


i 


\z — a J 


-1 


—i 



(6-1) 
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The difference from ^'o (5-1) is simply in the power growth near the endpoints a, a of the main arc. The 
transformation that links ^Pq to '^i is called (discrete) Schlesinger isomonodromic transformation. 

In fact the two matrices are simply related one to another as seen below 



1 


' -i -1 ■ 




i-A')a3 


' i 1 " 


2 


1 i 


\z — a J 




-1 -i 



Lemma 6.1 (Schlesinger chain) The matrices 

are related by a left-multiplication by a rational matrix 

^k{z) = i?K(z)*o(^), 

where 



K ez, 



Rk{z) 



2 



p 



(6-2) 



(6-3) 



(6-4) 



2 2 

Proof. The expression of Rk follows from straightforward computations; we only point out that the 
existence of such a rational left multiplier follows from the fact that all k solve the same RHP (jump 
conditions and normalization at infinity), but have different growth behaviors at the points a, a. Q.E.D. 

Mimicking the previous case of Definition 5.1, we shall state the following new definition. 

Definition 6.1 (Local parametrix near the spikes.) Let C,{z;£) be the local conformal coordinate in 
D, introduced in Definition 3.1, so that 

4 



-/i(z;x,t) =0(C;t) = -O+tQ'^ . 
e 5 

Let '^{^;v) denote the modified Psi-function (4-8) of the Painleve I problem with /3q 
/3_2 = ^1; 1^2 = ^^)- Then we define the parametrix 



Pl;a(z) 



1 

71 



1 -1 

i i 



T 3 



1 

-1 



.e(C;r)a3 



and we set 



V,,^{z) := {V,,„{z) f \ 



(6-5) 
= (and 

(6-6) 

(6-7) 



For brevity we will write simply Vi = Vi-a- We can then formulate the statement corresponding to 
Thm. 5.1 for the new local parametrix. 

Theorem 6.1 The matrix Vi satisfies: 

1. Withinlia, the matrix Vi{z) solves the exact jump conditions on the lenses and on the complemen- 
tary arc; 
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2. On the main arc (cut) Vi{z) satisfies 



Vi+{z) = cr27'i-(^)o-2 



(6-8) 



so that ^o'Pi within Dq, solves the exact jumps on all arcs contained therein ( the left-multiplier in 
the jump (6-8) cancels against the jump oJ"^q); 

3. The product "^o{z)'Pi{z) (and its inverse) are -as functions of z~ bounded within Dq, namely the 
matrix Vi{z) cancels the growth o/ at z — a; 



4- The restriction ofVi{z) on the boundary of Da is 



1 + 0(0 



(6-9) 



,1 + 

where ©(^"s) is uniform w.r.t. v in a neighborhood of a pole Vp not containing any zero of y{v). 

Proof. For the first three points the proceeds exactly as in Thm. 5.1 and hence is omitted. 

(4) Due to the asymptotic expansion (4-15) for ^I*, when restricted to the boundary 9Aq, we have 

(following the same computations as in Theorem 5.1) 

T 3 



^i(C) 



1 

1 



1 -1 

i i 

1 -1 
i i 



1 

-1 



.e(C)T3 



C 



-w (i + o(r^)) h 



i 1 

-i 1 

era 



(i + o(r^) 



VC + r/2-y 





t/2 + 




-C/y] 







(6-10) 



Q.E.D. 

Corresponding to this new local parametrix, we setup the approximation of the solution as 

£{z)'^i(z) for z outside of the disks 0^,0^, 

Y{z)= { f(z)*i(z)7'i;„(z) 
£iz)-^,iz)Vi,a{z) 



for z inside of the disk D^, 
for z inside of the disk 



(6-11) 



Recall that we are considering the regime v—Vp = 0{er,+'') (and hence y = 0{e- 



^)) where v 



3^2 



and v > Q. The boundaries of both disks Dc,,Dg are mapped by the conformal changes of coordinates 
C,C on some closed curves in the respective planes, that expand homothetically with a scale factor £~5. 
Then the behavior of the jump in (6-12) is determined by the behavior of the local parametrix Vi — Vi-a 
on the boundary of A^. The jump of £ is 



£+{z)=£^{z)^,iz)Vzi^^\z) , ze 



(6-12) 
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From eq. (6-9) it is clear that the rate of approach of Vi to the identity on the boundary is seriously 
impeded by the last factor in (6-9), which fails to converge to 1 when v — 0^ namely, when y — 0{e~i)). 
More precisely we have from (6-12) (since C,/y = ©(e^")) 

£+(z) =£_(z) (l + 0(£"""(s^2.)^ _ (g_^3) 

Before tackling the general problem = 0, we shall see what happens when v = Vp (namely y = oo) and 
we are exactly at the "top of a spike" 

6.1 The top of the spike: amplitude 

When V — Vp]s exactly a pole of the tritronquee we can use the expansion (4-16) for the expansion of the 
local parametrix on the boundary of the disks , . Since the first term after the identity is of order 
0{C^^) — 0{e^) when restricted on the boundary, the error term in (6-11) is then of the form 

f(z) = l + C'(£5,z-i). (6-14) 

Near z = oo we can write 0{ei , z^^) as and so we have 

Y{z) = E{z)-^,{z) ={l + -0{e^^ {^ + Yz 

We thus have 

q{x,t,e)=-2 lim zTi^2{z)^-2e^3(oo;x,t) z (£^r^) - ei*(^'*) f-36(a;, t) -I- ©(e^ )) =(6-16) 

=e;(*(^'*)-^-)36(l + 0(e5)), 
where {x, t) — v^^{vp) corresponds to the top of the spike. Since the diffeomorphism v = v{x, t, e) is scaled 

4 

as £5 , see Corollary 3.2, the corresponding to Vp spike will approach the point of gradient catastrophe 
{xo,to) at 0{e^) rate. Here b can be taken as the value at (xo,to) because the difference between the 
values at {xp,tp) and (a^oi^o) is of order 0{ei) (see (3-43) or (3-40)). Thus, we have proved the following 
theorem; 

Theorem 6.2 The asymptotic amplitude of a spike near the gradient-catastrophe point (xo,to) is (up to 
0{e^) accuracy) three times the amplitude predicted by the Whitham modulation equations at {xo,tQ). 

This result is a bit unexpected in that it entails a very simple universality; the three-fold amplitude 
of the first spikes appears to be entirely independent of the initial data. 

In fact the factor of 3 (together with the phase-shift, i.e., the minus sign) can be traced to the exponent 
— |cr3 of the outer parametrix 5*1, compared with the exponent of 5*0. In turn, this exponent is 
determined by the asymptotic behavior of the modified ^P-function for the PI problem (4-16). The latter 
is due to the shearing transformation for the cubic potential V{£,;a) appearing in (4-13). 



-0{z-'] 



(6-15) 
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It seems clear that, were we to study a nongeneric gradient catastrophe (i.e., more than one new 
main arc emerging from the endpoint of an already existing main arc), we would have to replace the 
PI problem by a higher member of the Painleve I hierarchy, which are characterized by exponents |, | 
etc. Thus, we conjecture that the amplitude of the first spikes after a nongeneric gradient catastrophe 
will be 5,7, 9, . . . times the amplitude at the gradient catastrophe {xo,to), depending on the degree of 
degeneracy. 

Conjecture 6.1 The amplitude of the spikes in the neighborhoods of the points of gradient catastrophe 
are odd multiples of the amplitude at the point itself. 

We reserve to verify this in a subsequent publication. The possible shape of these spikes in the case of a 
degenerate gradient catastrophe are discussed in Remark 6.1 below. 

6.2 The shape of the spike 

Due to the rightmost term in (6-9), the rate of approach of Vi to the identity on the boundary of Ci^a) 
in the regime 

V ~ Vp ^ 0{e^+'') ^y = 0{e-7^-'^'') (6-17) 

becomes slower as v approaches the critical value v — a.i which point the parametrix Vi does not tend 
to the identity at all. Indeed, since C — 0{£^i), we see that on the boundary of C(Aq) 




(6-18) 



so that the jumps of £ on the circles Aa, Ag are 

£+{z) = £_(z) (l + ©(e-^'^lMfej^s))^ . (6-19) 

From the standard approximation theorems for Riemann-Hilbert problems, it is seen that £ converges 
(uniformly) to the identity only up to the same rate of convergence of the jmnps, in this case ©(^"""(^(t^^s)). 
in particular the "error" becomes worse and worse as u ^ Q. 

As we shall presently see, it is possible (and necessary) to modify the outer parametrix in such a 
way that the troublesome factor Q{z) above is exactly taken care of. To account for the term Q{z) we 
shall seek an exact solution of the Riemann-Hilbert problem described hereafter: let C,{z)t(,*{z) — C(z) 
be the local conformal scaling coordinates in the neighborhoods of a, a respectively of the form 

C(z) =e"5C(z-a)(l + 0(z-a)) , l{z) ^ £-iC{z - a){l + 0{z - a), (6-20) 

2 

where C = (|Ci)^ 7^ 0. Let us assume that the circles 9Aq and 9Aa (oriented counterclockwise) have 
some small radius r, so that \C,{z)/y\ < 1 and \(*{z)/y\ < 1 respectively on these two circles. 
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Problem 6.1 Find a piecewise analytic matrix E{z) on the complement of the two circles described above 
and such that 



E{z) = l + 0{z-^) as z^oo, 

E+{z) = E^{z)^iiz)Q-\z)^^\z) = E^FM„{z)F~^ , \z - a\ 
E+{z) = E_{z)^i{z)Q-\z)^^\z) ^ E^FM^{z)F-^ , \z-a\ 

M^{z) ^ 
M^{z) 



v/1 - C/y 



j/(z-a)3 
1 



1 



y(z-a)^ 
1 



F 



1 

71 



-i -1 
1 i 



(6-21) 
(6-22) 
(6-23) 

(6-24) 

(6-25) 
(6-26) 



Before solving Problem 6.1 we show how its solution will be used. If E{z) solves Problem 6.1 then we 
re-define the approximation of Y as 



Y{z) 



for z outside of the disks Aq, Aq, 



£{z)E{z)-^i{z) 
£{z)E{z)-^ i{z)'Pi{z) for z inside of the disks Aq, A 
Then, according to (6-18), the jump of £{z) will be 



(6-27) 



E^E^^iViQ'if^'^EZ^ 



E+ = E-E^^iVi'H^^E^^ 

= £^E_ (l + Oiei)) EZ'=£-{1 + 0(e*) ) , 



(6-28) 

where the last equality holds provided that i?_(z), EZ^{z) are bounded uniformly in e on the bound- 
aries (which will be the case indeed). 

6.2.1 Solution to Problem 6.1 

The problem has an explicit solution. For the sake of simpler computations we will conjugate E{z) by 
the constant matrix F, so that E{z) := F^^E{z)F has the matrices Af^, M-^ for jumps. 
It is known that the solution E{z) (if it exists) must satisfy the integral equation 



E{z) 



E_{s){Mjs) - 1) ds 



\s—a\—r 



2iTT 



E_{s)iAUs)-l) ds 



\s—a\—r 



2i7r' 



(6-29) 



where E_{s) is the (analytic continuation of the) solution from the outside of the circles A^, Ag. It is 
crucial that the jump matrices Ma^a admit a simple decomposition of the form 



M„(z)-l = 0„(z) + 



1 




=: O^iz) + 
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-le 





1 



(6-30) 



where Oa (z) denote some locally analytic matrices in the respective neighborhoods (whose expression 
the reader can evince from the above formulae but which has no bearings in the considerations to follow) . 
What is essential in the following is that when evaluated at 2 = a and z ~ a, these local analytic matrices 
will be multiples of (t+ and (t_ respectively. 

Consider the Ansatz ^ 

^-(z)-l+^- + ^. (6-31) 

z — a z — a 

The expression of E{z) in the inside of the disks Aa,Ag (i.e. i?+(z)) has no particular interest for us 
and can be simply obtained from the jump condition. According to (6-31), the integral formula (6-29) 
becomes 



A <7^n„ 



z — a z — a 



We thus have the system 



A 



{a — z) {z — a)^ a — z {a — z){a — a) 
{a — z) (z — a)^ a — z {a — z){a — a) 



A(j+ = 0, 



A(T_ = 0, 



(6-32) 

(6-33) 
(6-34) 



The solution is given by 



A = 



(a— a)^ 



(a— a)^ 



Tla 



a — a 



2b\C\ 



1 



-i^ 




{2ibfC 
.|C|462 



2b\C\ 





£5 J, 



(6-35) 



Thus the solution of the problem for E has the form (6-31), (6-35) in the region outside the disks. 

6.2.2 Partial Schlesinger transformation and improved leading order asymptotics: shape 
of the spike 

Theorem 6.3 (Shape of the spikes) The spikes around the point of gradient catastrophe for q{x,t,s) 
are in one-to-one correspondence with the poles of the tritronquee solution y{v) of the Painleve I equation 



6y ~ V 



(6-36) 



47 



and have a universal shape. For a given pole v — Vp of the tritronquee solution, the shape of the spike is 
described by 



q{x,t,e) l + {s + 2i){s + 2i) 



qoixo,to,e) 



1+ s 



l + 0(£5)) 



- 3 + 4i3?(s) 
1 + kP 



1+0(65 



(6-37) 



where qo{x,t,e) is the genus zero approximation (see (2-20), (2-21), and Remark 3.7) the variable s is 
defined in terms of the tritronquee solution y{v) as 



5Ci ^ = 2ib 

2 



5Ci \ 5 2ih 



[v-vpf {l + 0{v-vpf) 



(6-38) 



and the diffeomorphism v — v{x, t, e) is given in Corollary 3.2. The formula and the error term are valid 
uniformly for (x,t) in a 0{£) -neighborhood of the center of the spike {xp,tp) = v~^{vp), or -which is the 
same- for v{x, t, e) — Vp = 0{ei). 

Proof. It is apparent that both E-{z) and EZ'^iz) are uniformly bounded on the boundaries of the disks 
Aq, Aa as long as y — 0{e~i~'^) for some u > Q. Thus, according to (6-28), the matrix £{z) will be 
1 -\- 0{£'s). Then the leading order approximation is thus given by 



*i = S(z)*i(z) = FE{z) 



z — a 
z — a 



where 



F 



1 

71 



-i -1 
1 i 



(6-39) 
(6-40) 



From this expression we find the behavior of the amplitude 

/ ^ \ 3 _ 

A A 



^i{z)^F \ 1 



z — a z — a \ z — a 



z — a 



3ib 

'2z 



-i 

1 



1 1 

1 + ^z 

£5 |y| 



F 



F-^ = 1 + F 



* — 2 



-ie-^J8iCh^/y 



A + A 3(a — a) 

z A z 

-ie--^^-8iCb^/y 

~* — 2 



Since q — — 21im^_>.tx) zTi2{z), where 

r(z) = ei's(°°)'"^£(z)*i(z)e-i's(°o)'^3^ 

we have 

(7(x,t,e) =e^5(°°) (-36 + 2(i^(A + A)^^-i)i2+0(e5)) = 



F-^ -f 0(45-^1) 



(6-42) 



(6-43) 



-36 



£i\y\ + 2b\C\ 



(|C|462 + iein (V-8jC63y)) + 0{ei)^ = (6-44) 



1 + 



2b\C\ 



e?*& 3 



2b\C\ 



4i5R 



-2ibC 



(6-45) 



48 



Since the variable w is a difFeomorphism on D, which is a size es neighborhood of the point of gradient 
catastrophe, then the variable s defined in (6-38) is a diffeomorphism on a neighborhood of size 0{e) 
around the spike. In terms of s we can rewrite (6-45) as 

q{x,t,e) = -e^ ^^-^^ ^(1 + O(£0)=e^ b ^-^^ (l + 0(£ = )). (6-46) 

Q.E.D. 

We remark in the following corollary that from 6-37 it follows immediately that each spike has two 
zeroes (nodes), one on each side. 

Corollary 6.1 For each spike near the gradient catastrophe point there are two zeroes (nodes) and they 
occur at the time t — tp (asymptotically as e ^ 0) and with x — Xp = ±e^\/3- 

6.2.3 The spike as a rational breather 

To gain a bit more insight into the formula (6-37) we start with the following observation; in terms of 
the coordinates t we have 

= -K- + 2{a + a)^V (l + 0{ei)) , (6-47) 

s = ±2ibe-^ {x-Xp + (4a + 2ib){t - tp)) (l + 0{e^)^ . (6-48) 

The phase can be expanded near the spike using (3-69) and retaining just the linear term; indeed the 
scale of the spike is 0(e) and hence x — Xp = 0{e) and t — tp = 0{e). We thus find 



So, 



where 



s = 2ib\^^ + 4ar] + 2ibr] j = -AVr^ + 2ib{£^ + 40??), (6-49) 

(6-50) 



j{x,t,e) = ei*(--*'')gfc,(^,77)(l + ©(e^)), (6-51) 



1 + Aib'^r] 



(6-52) 



462(^ + 4a?7)2 + m^rf ^ 

The expression Qbr{£,,il) is well known in the literature [24]; it is called the rational breather 
solution for NLS. Indeed it solves exactly the NLS equation 

tedtQbr + e^dlQbr + 2 \Qbr\^ Qbr - 0, (6-53) 

where 

?=^^, V='-^. (6-54) 
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In our case it is obtained from the "stationary" breather (depicted in Fig. 10) 

by applying the transformations (mapping sohitions into solutions) 

Qi^,Tj)^XQ{X^,X'ij), 0(^,77) =e^('=^-'='*)Q(^-2H,t). (6-56) 

The breather has the maximum equal to 3 at ^ = = ry and tends to 1 as \s\ — > oo; therefore it 
smoothly interpolates the regime near the spike with the regime "far" from it. 




Figure 10: The theoretical shape of the spike as a function oi s = i. 5- as given by Eq. (6-46). This 

V 

also represents the exact single-breather solution Qbr where s — 2ib{^ + Aarj + 2ibri). 

Remark 6.1 While completing the preparation of this paper we came across [1], where a hierarchy 0/ 
rational breathers generalizing the standard one (6-55) is investigated. Interestingly, the maximum am- 
plitude of these breathers is also an odd integer and hence this prompts the speculation that they should 
play the same role of (6-55) in the cases of gradient catastrophes with higher degeneracy (compare with 
Conjecture 6.1). 
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-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 



Figure 11: Some simple numerics (using Fast Fourier Transform) for the initial data 17(3;, 0) = 

and e = note that the amplitude at the first spike (numerically determined) is |(jpeafc| = 3.9256 and 
\qo\ = 1.3137. This is almost exactly three times, 3|go| — 3.9411. In this case ^ = and to — j as per 
Example 3.1. 



Remark 6.2 (Consistency Check) Suppose y remains fixed and e ^ in the formula; for A, A (6- 
35); we then obtain the limit 

(6-57) 



g — g 

Z — Q. 







F 



Z — OL 

s — g 











g — g 

Z — OL 



g — g 

Z — OL 



a 



z — a 



z ~ a 

-1-3 



F^i = F 



F-^=^o{z). 



(6-58) 



Example 6.1 According to numerical studies by [18], the closest pole to the origin of the tritronquee is 
at distance of about \vp\ ^ 2.38 (on the bisecant of the sector containing the poles). According to our 
Theorem 6.3, this is the first spike after the gradient catastrophe and its location is thus determined by 
the diffeomorphism v — v{x,t,e). 

In the case of the initial data of Example 3.1 with fi — 0, the solution is symmetric (x —x), hence 
the first spike occurs at x — xq ~ 0; the time of the first spike is then estimated at 



^ spike 



^to+l 



5|gi| 
4 



(26o) 



^|«,|(l + 0(et)) 



(6-59) 
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For the case of the initial data of Example 3.1 (fi — 0), the exact values of bo — b{xo,to) = \/2, ii 
and \Ci\ = were found in [26]; this gives 



1 




(6-60) 



In Fig. 11 we report on a simple numerical example with initial data 5(2;, 0) = sech{x), e = plotted 
are the profiles of\q{x,t)\ fort — to — l/A which is the theoretical time of gradient catastrophe; the second 
curve is the time where numerically the maximum amplitude is achieved at x — 0. The numerics indicates 
tspike — 0.28 whereas the formula (6-60) predicts tspike = 0.2791260482, which -considering that e is not 
that .small- is an amazingly close estimate. 

Our numerics was produced by integrating directly the NLS equation using FFT in the x-direction; it 
should be noted that the initial data we have chosen correspond to a pure soliton situation, and hence one 
could use a better numerical approach based on exact linear algebra as explained in [21]. 

A Estimate of the parametrix on a circle of large radius uni- 
formly for large y 

This section is rather technical; its goal is to prove that the sectorial analytic matrix function 
that is related to the solution P{£„v) of the RHP (4.1) through $(C;w) = G(f , t;)P(^, ^;)e''(^'")'^^ has 
asymptotic behavior (4-15) as ^ — >■ 00 uniformly with respect to v—Vp — 0(£5 ), where Vp is one of the poles 
of the tritronqucc solution y(v) of PI, or, equivalently, uniformly with respect to || {x, t)'~{xp, tp) \\— 0{e). 

In fact, we will prove a somewhat stronger result, stated in Theorem A.l below. Corollary A.l of this 
theorem has been used in the main body of the paper. 

Theorem A.l Let S be one of the sectors shown on Figure 8 and 5'(^, u) — G{S,,v)P{S,,v)e^'^^'^^'^'-^ be a 
solution of the ODE (4-9), obtained from the sectorial solutions P(^,w) of the Problem 4.I in the sector 
S. Let A be a constant nonsingular diagonal matrix and Vy be a disk of fixed radius r > centered at 
^ = y. Then, there exist some constants > 0, p2 > 0, both independent of y and of each other, such 
that \I'(^,w) has the representation 



1 




(A.l) 



^/2 



where 




A2 



(A.2) 



(A.3) 
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uniformly for S_ G S{(,^) = Sr\{£, : |^| > S,*}\Dy and for all sufficiently large y. Here = ^vy/^ 

and the contour of integration in (A. 2) is a ray in the direction of the bisector of S . In the case when 
y ^ 5^, the O term in (A.l) is independent of j^, i.e., it is 0(^^3). The same O term is valid even 
if y ^ S^_^ everywhere in S^^ except a certain "shadow" ofDy region, see Fig. 12 and description of this 
region below. 

Proof of this theorem constitutes the bulk of this Appendix. Let us consider sector S that is not 
adjacent to M_. This sector is bisected by the Stokes' ray (. where 3?^ 5 =0. Let us extend S by some angle 
5 € (0, |7r) in both directions. For definiteness, we shall focus on the sector S ~ {arg(^) e (—(5, ^ + (5)}, 
the remaining being treated similarly. It is more convenient to work with variable instead of ^; 

therefore, we introduce a new variable 



(A.4) 



We will consider u(^) as an analytic function (with the principal determinations of the roots) on sectors 
of the ^-plane of width up to 47r/5 and for |^| sufficiently large: therefore in any such region the map u{£) 
is holomorphically invertible because v will be chosen uniformly bounded (in fact- in a neighborhood of 
Vp), and we will denote its inverse by ^{u). We will write for brevity u or ^ with the understanding that 
they may be viewed as functions of the other variable. A direct calculation shows that 

2/5 

' (A.5) 



5i 



52/521/5 ^ ' 



We will let v — Vp he bounded above, say, by 2, and, since the branch-points of u(^) are at ±-y/-(;/2, the 
map u(^) is invertible in any sectorial domain {|^| > l^/v^ + 1|, as long as the opening of this domain is 
less than 27r. 

The image of S under the map (A.4) will be denoted by 6. Then, asymptotically for large u, 

& = {-(/) < arg(u) <(/)}, ^TT < < TT. (A. 6) 

It will be convenient for us to consider (A. 6) as the definition for (3, whereas S is the preimage of ©, i.e., 
S ~ u^^{&). By we denote sector 6 shifted from the origin to some uq G C, i.e., &uo = & + uo; for 
convenience we shall take uq E Without any loss of generality we assume uq to be sufficiently large. 
Our first step in proving Theorem A.l is to "prenormalize" the ODE (4-9), satisfied by ^, written with 
respect to the independent variable u. 

Lemma A.l Let W{u) be defined by 

4'(0 = 2'^3/2(4^3 _ 2^„)-3/4_L(^^ + ^y^mOWiO (A.7) 
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where the matrix T is defined by (A. 3). Then W satisfies the ODE in the variable u 



Wu = 



a-i + B{u) 



W, 



(A., 



where B{u) = 0{u^^^^) uniformly with respect to sufficiently large values of y in u & Ruo,y — \ ^■ 
Here A = u{Vy). 



Proof. The transformation 



V2 



reduces the system (4-9) to 
where 



M 



((73 + iG2) 



V'i . 3 



X 



2H 



M = 2Hj + Ma + Ms 
The three terms in M have uniform bounds in S^g^y = u^^{Rug,y), where uq = u{£,o): 

• Since ^ ^ D we have IM3I < 

• According to (4-7) and the fact that H'{v) = y{v), we have 

2Hi = 2Hiv) - = 28/3 + 0(w - Vp); 

y(v) 



(A.9) 



(A. 10) 



(A.ll) 



(A.12) 



According to (4-7), (Iny)' — 2y^(1 -|- 0{v — Vp)'^) for ah sufficiently large y, so we now need to 

< ^ outside the domain \£,-y\ < dol^l- Thus, M2 = 0(^5) outside 



estimate Note that 
this domain. Inside the latter domain but outside Vy we have |4| = 0{£,^), so that 

M2 = 0(^5) 

for all ^ e S^g,y and uniformly in all sufhciently large y. 
Rewriting (A. 10) in the variable u we obtain 



(A.13) 



A",, 



.J2Vl_ 



M 



-Ox 



2(4^3-2^1;) 2 ^4^3 - 2iv 

6^2 _ M 



(0-3 + 10-2) 



X 



2(4^3 -2^w)3/2'^^ ' 2(4^3 _25i; 



(f73+«'T2)+0(u~"/'"') 



X. 



(A.14) 
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Dropping all the terms in M except M2, the previous equation becomes (see the buUeted list above) 

eC^-V _ , y'C(CT3+iCT2) 



Xu = i 
Using (A. 13) we have 



CT3 



(i + o(r')) 



X. 



0{u- 



and hence 



Xu = i 



2(4^3 - 2e^;)3/2^^ ' 8ey{^-y) 

for u G i?«o,y uniformly in y (for sufficiently large y). Setting X — TW with T as indicated in (A. 3, A. 7), 
simplifying and keeping track of the orders already estimated, yields the equation (A. 8). Q.E.D. 

We will use the ODE (A. 8) in order to setup a recursive approximation scheme within each sector; 
the key is that the estimate for B{u) in (A. 8) is uniform for w in a neighborhood of the pole Vp (hence, 
for large values of y) and |^ — y\ bounded below. 

The domain A — u{Dy) in the u plane is of size ^^([^(j/)!^). For simplicity (and without real loss of 
generality) we modify D to be the preimage of the exact disk A in the w-plane centered at u{y) of the 
radius do\u{y)\i . 



y (0-3 + icr2) 



X 



(A.15) 



(A.16) 



(A. 17) 



Lemma A. 2 Any solution W{u) to (A. 8) can he written in the form W = ^e'l'"^'^^ where 



du 



q(u) 



i + 



ly'i 



&eyii ~ y) 



and $(u) solves 



5'' = q'[cr3,*] 
with B = B{u) as in Lemma A.l (eq. (A. 8)). 

Proof. A direct substitution of the proposed expression into (A. 



(A.18) 
(A.19) 



Q.E.D. 



The differential equation (A. 8) is equivalent to the following Volterra integral equation 

$(u) = A + e'l^'') ( / e-i(")B$e'i('')d77 I g-'^")'"^ = A + ![$](«), 
\Jn(u) J 



where q from (A.18) can be written as 

q(u) : 



ly 



du 



8y J 00 {^{u)~y)e{uy 



(A.20) 



(A.21) 



The symbol Jj^^^-j denotes the integration along a set of contours originating at u and extending to u = 00, 
with a different direction of contour Vli j{u) for different entries of the matrix, see Fig. 12 for the choice of 
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dO 








AO 


( uiyy\ 




Wo \ 




\A J 










- 










Figure 12: Region i?MQ,3^ is the complement of A and the union of A° (darker) and y (lighter). Shown 
in the picture are the contours of integration ilij{u), u G Ruo,y^ ihj) entry of the Volterra 

operator (A. 20). 



the contours (the exponential of the integrand should decrease in the direction of the contour, traversed 
from u to oo). It is promptly seen from the fundamental theorem of calculus is equivalent to (A. 18). The 
matrix A is a constant of integration, which, at this point, we choose to be a diagonal and nonsingular 
matrix, but otherwise undetermined. 

Changing variable u — u{£,) in the integral (A. 21) and using ^ ~ — 2z^^/^ + ^ and 

- Vln#^, (A.22) 



Viit-y) ^ VI- 
we find 

q(.) = - + ^ In J;'^). + 0{u-"') (A.23) 
in the region Ruo.y uniformly in y. Note that 

-4 = ^(x; - v,)^ + 0{{v - V,)') = -l+ 0{y-^) (A.24) 

and hence 

eq(uV3 =e(t«^-<*)-3 f ^^+^ Wi^o{u-^'\y-^)) (A.25) 

V Vi-y J 

in the region Rua,y uniformly for large y; the cut of the logarithm is taken from y to ^oo parallel to the 
Stokes line (bisecant to the sector). 
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The strategy is to show that the integral operator T, defined in (A. 20) is contractive in a suitable 
space, so that a solution can be sought through an iteration scheme; the analysis depends on the region 
u belongs to. For practical reasons we state the following lemma. 



Lemma A. 3 The 



function {^-^^^^^^ and its inverse satisfy the inequalities 



outside the circle y\> r, with Kq a constant independent of y (but depending on r). 



(A.26) 



Proof. The function ^^'^^^^^ is continuous (in fact smooth except at ^ = and ^ — y) on C\ {y} and 
has limit 1 at ^ = c», being unbounded in any neighborhood oi ^ = y. It can be also seen by calculus 
that it has no critical values in C \ {0, y} and hence the maximum is attained on a disk around y. The 
value there is easily estimated as in (A.26). Similar arguments work for the inverse function. Q.E.D. 

A.l Convergence of iterations in Ru^^y 

The region i?,^^ y is the lighter shaded region in Fig. 12; it is obtained by excising from Ruo,y the "cone 
of shadow" of the disk A with the angle ip2 G (0, tt) indicated in Fig. 12 chosen arbitrarily and fixed once 
and for all. The region A° will be such cone of shadow. For a given u S R^^ y, the contours of integration 
riij(u) are also shown in Fig. 12. It is important that the collection of contours fl{u) C R'^^^ y for any 

Lemma A. 4 // a matrix-function x(u) satisfies \\x{u)\\ < c|u|~™ in i?^^ y for some c > and some 
m > then \\Xix{u)\\ < cKi\u\^"^^^^^ , where Ii denoted the Volterra integral operator (A. 20) with the 
contours Q.{u) specified above. The constant Ki does not depend on x, uq G M+ and y G C, but depend 
on ni (it was assumed above that Uq and \y\ are large). 

Proof. The diagonal entries of Bx are unaffected by the conjugation by e'''^^ and are of order 
after integration they become of order u^™^5 . Regarding the off-diagonal entries of Bx, let F{u) denote 
one of them; then, according to (A. 20), |-F('«)| < ci|u|~™^5 on Ruo.yj where the constants L > and 
TO > do not depend on y and we have (from (A. 25)) 

2 



.±2q(«; 



F{u) 



„±2iu 



F{u) 



iI + Oiu-'/\y-')), 



so that 



(/ + 0(.-3/^y-)) 



(A.27) 



(A.28) 



Due to the choice of the contours Q{u), the integration along the corresponding Qi j{u) can change the 
previous estimate only by a constant that does not depend on uq and y (but depends on 02)- Q.E.D. 
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Let us consider successive iterations <&„ = Ii$„_i, n > 1, of equation (A. 20), where $o = 0. Let 
A$„ = $„ — Then $i — A<i>i — A and A$„ = ZiA<i>„_i. We are now going to prove uniform in y 

convergence of the series in different subregions of Ruo,yi provided that uq, which is independent 

of y, is large enough. In all the analysis below y is assumed to be large. 

Let ||A|| = c. Then, according to Lemma A. 4, ||A$2|| < cKi\u\~'^/^ in y uniformly in uq and in 
y. We can choose uq G M+ so large that Ki\u\^'^/^ < \ for all u e Ruo,y This choice of uq guarantees 
convergence of the series J2T in Ruo,y uniformly in y. Moreover, we obtain the estimate 

= A[/ + 0(u-i/5)] (A 29) 

in R'^^ y uniformly in y. 

Let Ti-^o.y denote the preimage of y under the map u = it(^), where uq — u{£^q). Going back to 
the system (4-9), we obtain the following result. 

Lemma A. 5 The solution w) to the system (4-9) has behavior 

4f(^, v) = 2'^^/2(4^3 _ 2^i;)-'^3/4 1 ^ (T3)rA(/ + 0(r5))ei("(«» (A.30) 

v2 

in TZ^g,y uniformly in y provided is large enough. Moreover, for a fixed diagonal A, condition (A.30) 
uniquely determines the solution ^(^,?;) of (4-9). 

Comparison of (A.30) and (4-8) together with the asymptotics of "if (4-6) and (4-14) yields 

A = diag(l,i). (A.31) 

A. 2 Convergence of iterations in A*^ 

If y does not belong to the sector S we are considering, we can always arrange that Vy n S* = 0, so that 
^5o,y coincides with the preimage of &uo and, thus, the estimate of Lemma A. 5 holds throughout that 
sector. We want to extend the statement of Lemma A. 5 into the preimage 2?° of the region A° for the 
case that y £ S. Because of the construction of S we can assume, without any loss of generality, that 
y G S, where S is any proper subsector of S. 

Consider the tangent line A to the disk A (which is centered at u{y)) that is parallel to e""^"^ and 
located above A, i.e., if is the point of tangency, then argit* > argu(j/). Here £ ((j),(f)2)- It is clear 
that A divides &uo into two regions. Let i?^^ y denote one of these regions, namely, the one that does 
not contain A. We want to extend the statement of Lemma A. 5 into the region i?^^^ y. To this end, we 
construct a solution $ of the integral equation (A. 20) in i?^^^ y by successive iterations. Let ui denote 
the point of intersection of 96„„ and A. The collection of contours of integration il{u) in the case of 
the region i?^ is similar to the case of the region i?^ , considered above, except that the contour of 
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integration for the entry (1,2) is the segment [ui,u], provided that [ui,u] C Ruf,,y For u € Rug^y that 
do not satisfy the latter condition, the contour of integration is the union of [wi,uo] U [uo,u]. 

Let I2 denote the integral operator in the Volterra equation (A. 20) in the region i?^^ ^ with the 
contours defined above. We can repeat the previous estimates of integrals to extend the statement of 
Lemma A. 4 to the integral operator I2 in the region i?^^ y-. the only difference is the finite contour of 
integration, where the desired estimate comes from Lemma 14.2, [29]. The solution $ to (A. 20) thus 
satisfies the same estimate of Lemma A. 5. We also have 

l>(u) = A[/ + 0(u-i/5)] (A 32) 

in R!^^ y uniformly in y. Here A is given by (A. 31). 

Let mq G M_|_ be fixed: for any sufficiently large \u{y)\ the set y n y consists of two disjoint 
components y shown in Fig. 13 (otherwise, we can extend A in such a way that it will intersect the 
boundary of ©uq ^-nd the corresponding ISP disappear). Let S{y) denote the Stokes matrix connecting 
solutions of (A. 8) ^e'^'^^ and the previously constructed one ^e"^"^ . Then 

cI>(^i,y) = |.(u,y)e'^("■^)'^^5(J/)e-'^("■^)-^ ^(y) = J '^^^ ^ (A.33) 

where the triangularity of S{y) follows from the fact that A as m cx), u G T^^^ y. Writing 

q = qiu^y) in (A.33), we emphasize that q depends on y. 

Estimate for s(y). We now want to obtain an estimate for s{y) by comparing $ and $ in the region 
TJj^ y. Let us denote by $ the analytic continuation of $ from T^^^ ^ throughout R^^^y (throughout the 
dark shaded region) into T^^^ y. According to (4-8), '^{£,,y) has monodromy —1 as ^ goes around C ~ V- 
Because of (A. 7), W = $e'i'^^ has the same monodromy. It follows from (A. 23) that e'"'^^ has monodromy 

Mniy :— exp 2i7r-^(T3 around ^ = y. Therefore we have 

4y2 



$(u; y) = -$(u; y)m-' , u G (A.34) 



The matrix niy (which is constant in ul) has the behavior niy = 1 + 0{y ^) on account of (A. 24). 
According to (A. 29 and the fact that niy is uniformly bounded (for all large y), we have 

ll$-A|| <i^2|^^r^ \\^-A\\<K2\u\--^ iiiTl^^y, (A.35) 

with a constant K2 independent of y. But, according to (A. 32), $ too satisfies the same estimate. 
Substituting these estimates into (A.33), we obtain 

s(j;)e2i(") =0(^-5) (A.36) 
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Figure 13: The region Ru^^y (shaded) and its two subregions T^^ ^ (hghter shading). The goal is 

to extend the estimate (A. 30) into the remaining third (darker shading) subregion of Ru^^y 



in Tjj^ y uniformly in y. In particular, according to (A. 23), \s{y)\ < i^ajuilse^^^^ for some > that 
is independent of y. Therefore, in the region A° U R^^^ y (darker shade on Fig. 13), we have 



s{y)e 



2q(«) 



2Q{ui—u) 



(A.37) 



Given the geometry of regions and the fact that y G S, we have 

9(u - Ml) > Kiluiy) -ui\> K5\u{y) - uo\, Vu G A°, 



(A.38) 



where K^, denote some positive constants independent of uq, y. Thus, the matrix ei("'^)'^='S'(y)e~'''^"'^^'^^ 
in (A. 33) is exponentially close to 1 as y — ^ oo, where u{^q) = uq. So, according to (A. 33), the 
estimate of Lemma A. 5 (with a modified O term) can be extended to A° (more precisely, to the shaded 
part of it A° n Rug,y), provided y — ^ oo. Because of the construction of 5", there exists /i > 0, such 
that S contains S'(/i|'Col)i where by S now we mean the original sector of the RHP (8), see Theorem (A.l). 
We denote ^* — fJ,\£,o\. Existence of the positive constant p2 in the estimate (A.l) follows from (A. 33) 
(A.37) and (A.38). The case of sectors S adjacent to the jump contour M_ on Fig. 8 can be considered 
similarly since we can rotate this contour in both directions by angle up to ^. 
The proof of Theorem A.l is completed. 

According to (A. 3) and (A. 13), the matrix T from (A.l) can be absorbed into 1 + 0(^~3) term. To 
within the same estimate we can replace also \/4^^^~^-^ by 2^5 and recast the theorem into 
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Corollary A.l Under the same assumptions and in the same notations of Theorem A.l we have 

*(e,^')=^^-^(-l+a3)(A + 0(^^y-^e--l*l^^^))(^^^)"e-(«-)- (A.39) 
Theorem A.l together with the convergence of iterations imply the following corollary. 
Corollary A. 2 For any sufficiently large ^ we have Imiy^y^ — 'I'(^,Up). 
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